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ABSTRACT 


A  computer  program  (written  in  FORTRAN)  is  presented 
which  uses  Monte  Carlo  techniques  to  simulate  one-searcher, 
one-target  passive  acoustic  ASW  search  that  terminates  at 
detection.  A  threshold  crossing  detection  model  is  used, 
and  stochastic  variations  in  the  acoustic  signal  are  modeled 
using  either  a  Lambda- Sigma  Jump  or  Gauss-Markov  error 
process.  Both  platforms  have  the  capability  of  detecting 
each  other,  and  area  and  barrier  searches  are  modeled. 
Features  of  the  program  include  interactive  data  input, 
extensive  use  of  graphical  displays,  and  thorough  statis¬ 
tical  analysis  of  the  results  of  the  simulation. 
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I.  DESCRIPTION  OF  THE  SIMULATION  MODEL 

A.  INTRODUCTION 

Although  the  body  of  analysis  devoted  to  search  theory 
is  extensive  (see,  for  instance,  the  lists  of  references  in 
Koopman  [Ref.  1]  and  Washburn  [Ref.  2]),  most  of  the  tech¬ 
niques  are  applicable  only  to  stationary  targets.  Recently 
there  have  been  some  important  advances  in  the  theory  of 
search  for  moving  targets  (see  Eagle  [Ref.  3],  Brown  (Ref.  4], 
Stewart  [Ref.  5],  Stone  [Ref.  6],  and  Washburn  [Ref.  7]). 

Even  so,  it  is  still  a  very  difficult,  if  not  impossible, 
computational  problem  to  calculate  the  probability  of 
detection  or  the  mean  time  to  detection  when  the  searcher 
has  a  speed  close  to  that  of  the  target  and  follows  a 
realistic  (i.e.  not  random)  track.  As  a  result,  computer 
simulation  of  the  search  process  is  often  used  to  evaluate 
the  effectiveness  of  search  tactics. 

The  Monte  Carlo  simulation  computer  program  presented  in 
this  thesis  is  an  attempt  to  model  the  passive  acoustic 
search  process  in  a  one-searcher,  one-target  scenario  where 
both  the  searcher  and  target  platforms  have  sensors  capable 
of  detecting  each  other.  The  results  of  the  simulation 
allows  the  analyst  to  closely  approximate  parameters  (e.g. 
cumulative  probability  of  detection,  detection  times,  detec¬ 
tion  ranges)  which  are  essential  in  planning  search 

operations.  With  these  results,  the  military  planner  can 
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make  informed  decisions  concerning  asset  allocation  and 
specific  platform  employment. 

The  computer  program,  Passive  Acoustic  Search  Simulation 
(PASS) ,  is  based  on  an  acoustic  fluctuation  model  described 
by  McCabe  [Ref.  10:pp.7-23]  and  Hurley  [Ref.  lltpp. 25-34] . 
Precedessor  programs  include  the  FORTRAN  program  BEAR,  and 
the  ALGOL  program  SEARCH.  BEAR  was  written  by  L.  K.  Arnold 
of  D.  H.  Wagner,  Associates  (see  Arnold  [Ref.  8]  for  program 
documentation) .  SEARCH,  a  derivitive  of  BEAR,  was  written 
by  W.  J.  Browning  and  J.  Risberg,  also  of  Wagner, 

Associates  (see  Browning  [Ref.  9]  for  program  documentation) . 
SEARCH  has  been  used  extensively  in  operational  and  exercise 
environments  under  the  direction  of  Submarine  Development 
Squadron  Twelve  (SUBDEVRON-12) ,  where  it  is  currently 
resident  on  the  Semi-Automatic  Reconstruction  Facility 
(SARF)  computer.  The  basic  structure  of  PASS  is  similar  to 
that  of  BEAR,  with  some  of  the  features  of  SEARCH  added. 
Additionally,  some  modifications  were  made  to  algorithms  in 
these  earlier  programs,  and  new  features  were  added, 
including  interactive  data  input  and  extensive  use  of 
graphical  presentation  of  simulation  results. 

PASS  is  written  in  FORTRAN  and  is  designed  for  use  at 
the  Naval  Postgraduate  School  (NPGS) .  PASS  is  not  portable 
in  that  it  uses  external  subroutines  in  the  Non- International 
Mathematics  and  Statistics  Library  (NONIMSL)  and  the  DISSPLA 
graphics  system.  The  -  subsequent  sections  of  this  chapter 
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are  designed  to  give  the  reader  an  overview  of  the  model 
without  going  into  the  specifics  of  the  program  code  or 
logic.  If  the  reader  is  interested  only  in  learning  how  to 
use  PASS ,  it  is  recommended  that  this  chapter  and  Appendix  A 
be  studied.  Appendix  A  is  designed  as  a  User's  Guide  and 
does  not  require  any  knowledge  of  the  program  code  other 
than  variable  names,  which  are  described  in  Appendix  B.  If 
the  reader  is  interested  in  the  structure  of  the  program  and 
the  model,  then  this  chapter  and  Chapter  II  should  be  studied. 
Some  examples  of  varied  applications  of  PASS  are  shown  in 
Chapter  III. 

B.  GENERAL  PROGRAM  DESCRIPTION 

Program  PASS  is  a  Monte  Carlo  simulation  of  two  passive 
sonar  platforms  in  a  search  scenario.  Arbitrarily,  one 
platform  is  called  the  "searcher",  and  one  the  "target", 
where  the  target  is  typically  thought  of  as  a  submarine,  and 
the  searcher  as  any  passive  sonar  platform.  Since  both  ships 
have  passive  sensors,  the  assignment  of  the  names  searcher 
and  target  does  not  connotate  strict  roles  (i.e.,  the  searcher 
may  be  detected  first,  and  becomes  a  target  in  the  truest 
sense l ) . 

During  each  replication,  the  target's  motion  is 
partially  random  within  a  user-defined  search  area.  The 
searcher  motion  is  deterministic  and  controlled  by  user 
input.  Each  sensor  is  subject  to  random  acoustic  fluctuations 


which  induce  random  variations  in  detection  ranges  over 
time.  Detection  occurs  whenever  one  platform  obtains  a 
positive  signal  excess,  unless  the  user  employs  an  integra¬ 
tion  time  model,  which  requires  weak  signals  to  be  present 
longer  than  strong  signals  for  detection  to  occur. 

The  output  of  the  program  consists  of  distributions  of 
parameters  of  interest  (e.g.  time  of  detection  and  range  of 
detection) ,  and  point  estimates  derived  from  the  statistics 
of  these  distributions  (e.g.  mean  time  to  detection  and 
mean  detection  range) .  The  format  of  the  output  is  numerical 
for  point  estimates,  and  graphical  for  distributions. 

C.  SEARCH  SCENARIOS  USED  IN  PASS 

1.  Area  Search 

The  target  is  constrained  to  move  within  a  specific 
rectangular  area  ("search  area”)  in  a  "semi-random"  fashion. 
Target  speeds  are  distributed  uniformly  over  a  user  speci¬ 
fied  speed  range,  and  courses  are  uniform  on  (0,  360) 
degrees.  Times  between  course  and  speed  changes  are  exponen¬ 
tially  distributed.  The  size  of  the  search  area  and  the 
parameters  of  the  random  motion  are  user  specified.  The 
searcher  moves  along  a  track  which  need  not  be  contained 
within  the  search  area.  The  searcher  track  and  speeds  are 
user  specified. 

2 .  Barrier  Search 

The  target  is  constrained  laterally  to  be  within 
certain  user  specified  bounds,  thus  simulating  a  "choke  point" 


Target  motion  through  the  choke  point  is  generally  vertical 
(i.e.  from  "top"  to  "bottom")  with  random  course  and  speed 
deviations  parameterized  by  the  user.  The  searcher  track 
is  specified  as  in  Area  Search. 

The  distribution  of  the  initial  lateral  position  of 
the  target  can  be  specified  by  the  user  as  either  uniform 
across  the  choke  point,  or  with  specific  probabilities  as  a 
function  of  distance  across  the  choke  point.  In  using  the 
latter  option,  and  by  limiting  the  maximum  course  deviation, 
the  user  can  control  the  distribution  of  the  lateral  target 
position  when  the  target  penetrates  the  barrier. 

D.  THE  ENVIRONMENTAL  MODEL 

The  environment  is  modeled  in  a  conventional  manner. 

PASS  requires  range/propagation  loss  information  for  both 
the  searcher  and  target  sensors.  Direct  path  (DP)  and 
convergence  zone  (CZ)  data  are  entered  separately.  A  cubic 
interpolation  routine  is  used  to  determine  propagation  loss 
between  the  user- input  data  points.  The  convergence  zones 
are  modeled  as  inverted  "square-wells"  superimposed  over  the 
DP  propagation  loss  curve.  Up  to  five  convergence  zones  are 
allowed. 

E.  THE  DETECTION  MODEL 

In  PASS  a  detection  is  always  a  "secure  detection".  That 
is,  a  detection  by  the  searcher  means  that  the  searcher 
detected  the  target  before  the  target  detected  the  searcher. 
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Simultaneous  detections  are  treated  as  special  cases  where 
neither  the  searcher  or  the  target  gets  credit  for  a 
detection.  Hie  simulation  terminates  (replication  ends) 
whenever: 

1.  Hie  searcher  makes  a  detection. 

2.  Hie  target  makes  a  detection  (sometimes  referred  to 
as  a  counter-detection) . 

3.  A  simultaneous  detection  occurs. 

4.  In  the  barrier  scenario,  the  target  crosses  a  user 
specified  lower  boundary. 

5.  User  specified  maximum  search  time  is  exceeded. 

The  propagation  loss  curves  are  used  in  conjunction  with 
searcher  and  target  figure-of-merit  (FOM)  to  determine 
detection  ranges.  The  program  uses  FOM  data  which  must  be 
precalculated  by  the  user  based  on  specific  platform  source 
levels  (Ls) f  backgroung  noise  (self  noise  and  ambient  noise) 
and  sensor  directivity  (Le) ,  and  processor  (machine  and 
operator)  recognition  differential  (Nrd) . 

1.  Searcher  Figure-of-Merit 

The  Ls  of  the  target  is  assumed  to  be  independent 
of  target  speed.  In  effect  flow  related  noise  and  speed 
dependent  propulsion  noise  are  ignored.  At  the  user's 
option,  target  Ls  may  be  made  a  function  of  target  aspect. 

It  is  assumed  that  the  deterministic  sensor  performance  of 
the  searcher  is  affected  by  the  searcher  speed.  This  results 
in  a  searcher  FOM  that  is  dependent  on  searcher  speed,  and, 
optionally,  target  aspect. 
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2.  Target  Figure-of-Merit 

The  Ls  of  the  searcher  is  assumed  to  be  speed  depen¬ 
dent,  and  independent  of  searcher  aspect.  This  results  in  a 
target  FOM  that  is  dependent  on  searcher  speed. 

3.  The  Threshold  Crossing  Model  With  a  Stochastic  Error 

Process 

Given  a  FOM  against  a  target  for  a  search  platform, 
the  maximum  detection  range  of  the  searcher  is  determined 
from  the  propagation  loss  curves .  This  results  in  the  sensor 
conforming  to  the  "definite  range  law",  or  "cookie-cutter" 
performance.  This  type  of  sensor  readily  lends  itself  to 
mathematical  analysis,  but  has  the  disadvantage  of  not 
reflecting  actual  sensor  performance.  What  is  seen  under 
operational  conditions  is  a  variation  in  signal  strength 
(and  thus  a  variation  in  detection  ranges)  over  time,  even 
for  targets  at  fixed  range  and  source  level. 

To  model  these  variations  in  received  acoustic 
signal,  a  stochastic  error  process  is  added  to  the  determin¬ 
istic  figure-of-merit.  This  results  in  the  signal  excess 
consisting  of  deterministic  and  random  components.  For  a 
detailed  discussion  of  threshold  crossing  models  and 
associated  stochastic  error  processes,  see  McCabe  [Ref.  10] 
and  Hurley  [Ref.  11].  See  Appendix  D  for  a  discussion  of 
the  relationship  of  signal  excess  (SE)  to  figure-of-merit 
(  FOM)  . 

The  mean  signal  excess  (SE)  is  defined  as  the 
expected  difference  between  the  signal-to-noise  ratio  (SNR) 
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at  the  hydrophone  output  (in  decibels,  db)  and  the  SNR  that 
is  estimated  to  afford  a  probability  of  detection  of  0.5. 

The  model  used  for  the  "observed"  signal  excess,  as 
a  function  of  time  in  PASS,  is  of  the  form: 

SE  ■  31  +  X  +  Y  (1.1) 

where  SE  is  the  deterministic  mean  signal  excess,  and  X  and 
Y  are  random  variables  that  describe  the  fluctuations  in  the 
signal  level. 

If  we  let  the  subscript  s  apply  to  the  searcher  and 
t  for  target,  the  model  used  in  PASS  is: 

SEg(t)  -  SEg (t)  +  X(t)  +  Y(t)  (1.2) 

SEfc(t)  *  SEfc(t)  +  Y(t)  +  Z(t)  (1.3) 

where : 

1.  X ( t)  is  a  stochastic  error  process  which  describes  the 
fluctuations  "local"  to  the  searcher. 

2.  Z (t)  is  a  stochastic  error  process  which  describes  the 
fluctuations  "local"  to  the  target. 

3.  Y ( t)  is  a  stochastic  error  process  which  describes 
acoustic  fluctuations  of  a  global  nature  whose  effect 
is  common  to  the  searcher  and  target. 

4.  X ( t) ,  Y ( t) ,  and  Z(t)  are  independent  processes. 

X(t)  and  Z (t)  may  be  thought  of  as  modeling  the 
fluctuations  that  affect  only  the  performance  of  the 
detecting  platform,  such  as  onboard  fluctuations  in  sensor/ 
processor  performance,  including  operator  performance. 
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Y ( t)  may  be  thought  of  as  modeling  the  gross 
acoustic  fluctuations  that  occur  along  the  "environmental 
path”  between  the  platforms,  such  as  temperature  and 
salinity  gradients  and  changes  in  bottom  conditions. 

With  signal  excess  now  a  function  of  the  mean  signal 
excess  and  the  stochastic  effor  processes,  we  have  the  means 
to  model  the  behavior  of  the  sensor  by  employing  a  threshold 
crossing  detection  model  as  follows : 

1.  If  SE(t)  is  less  them  0,  the  probability  of  detection 
is  0.0. 

2.  If  SE (t)  is  greater  them  or  equal  to  0,  the  probability 
of  detection  is  1.0. 

Thus,  a  detection  occurs  at  time  t  if,  and  only  if,  SE(t)  is 
greater  than  or  equal  to  zero.  The  fact  that  SE(t)  is  now 
undergoing  random  fluctuations  about  the  mean  signal  excess 
results  in  random  variations  in  the  sensor  detection  range. 

F.  STOCHASTIC  ERROR  PROCESSES 

PASS  allows  the  user  to  choose  either  a  Lambda-Sigma 
Jump  (LSJ)  or  a  Gauss-Markov  (GMA)  error  process  to  model 
the  acoustic  fluctuations.  Hurley  [Ref.  11]  gives  an  over¬ 
view  of  these  processes,  and  others  used  in  search  simulation 
programs.  McCabe  [Ref.  10]  provides  a  more  mathematically 
complete  treatment  of  the  processes ,  along  with  a  compre¬ 
hensive  comparison  of  the  LSJ  and  GMA  models.  In  PASS  we 
use  the  sum  of  two  LSJ  or  GMA  process,  frequently  referred 

m 

to  as  a  compound  process.  In  the  following  discussion,  a 
simple  process,  vice  the  compound  process  is  considered  for 
clarity. 
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1.  The  Lambda- Sigma  Jump  Error  Process 

Define  a  Lambda-Sigma  Jump  error  process,  £(t),  as 

follows : 

5<t+s)  =  z(t)C(t)  +  [l-z(t)]n  (1.4a) 

(s^t) 

{1  if  s  <  T 

(1.4b) 

0  if  s  *  t 

where : 

1.  t  is  an  exponential  random  variable  with  rate  parameter  X. 

2.  n  is  a  normal  random  variable  with  zero  mean  and 
variance  a2 

Hie  process  can  be  though  of  as  a  "marked"  Poisson 
process  with  rate  parameter  X,  where  the  magnitude  of  the 
marks  are  themselves  random  variables  with  a  normal  distri¬ 
bution  (with  zero  mean  and  variance  a2) .  The  standard 
deviation  of  the  normal  distribution,  a,  can  be  thought  of 
as  a  scale  parameter  in  this  process.  The  figure-of-merit, 
then,  is  constant  over  exponentially  spaced  time  intervals, 
and  changes  to  new  levels  determined  by  a  normal  distribution 
at  the  end  of  each  of  these  time  intervals. 

The  covariance  function  for  £(t),  is: 

Cov  IC(t),  C(t+s) 3  =  o2e"Xs  (1.5) 

This  function  is  a  measure  of  the  correlation  between  values 
of  C(t)  at  different  times.  From  this  we  can  see  that  £(t) 
is  second  order  stationary.  That  is,  the  covariance  depends 
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on  the  time  difference,  s,  and  not  the  time,  t.  The  para¬ 
meter  X  determines  the  amount  of  dependence  between  successive 
acoustic  levels,  and  can  be  chosen  to  yield  any  result  from 
independence  to  complete  dependence.  The  resulting  process 
is  then  a  piecewise  continuous  function  in  which  there  are 
periods  of  complete  dependence  of  detection  opportunity 
interrupted  by  fluctuations  that  introduce  independence 
between  these  periods. 

In  PASS,  then,  we  must  provide  a  total  of  six  param¬ 
eters  to  completely  describe  the  compound  error  processes 
used  to  model  the  acoustic  fluctuations.  They  are: 

1.  X1=rate  parameter  for  the  searcher-local  error 
process,  X(t) . 

2.  ax=  scale  parameter  for  the  searcher- local  error 
process,  x(t)  . 

3.  X2*rate  parameter  for  the  global  error  process,  Y(t). 

4.  a2 = scale  parameter  for  the  global  error  process,  Y(t). 

5 .  X  3  *  rate  parameter  for  the  target  local  error 
process,  Z(t). 

6.  ct3=  scale  parameter  for  the  target  local  error 
process,  Z(t). 

The  selection  of  appropriate  values  for  X  and  a  is  a 
not  well  understood  function  of  the  environment,  sensor,  and 
processor.  In  practice,  the  selection  of  the  values  are 
subjective  decisions  based  on  the  experience  of  the  analysist 
as  what  seems  to  yield  reasonable  results.  Tehan  [Ref.  12] 
provides  a  discussion  of  estimation  techniques  for 
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these  parameters.  In  applying  these  recommendations  to 
PASS,  which  has  compound  processes,  the  values  of  X  and  a 
were  governed  by: 

1 

3^+X^S  6< (a  1+02)  2-9 

(1.6)  2  i  (1.7) 

3<X  +\  <5  6 <(a  +a  )  <9 

2  3  2  3 

See  Appendix  E  for  a  description  of  the  computer 

simulation  of  LSJ  McCabe  [Ref.  10:pp.9-13]  gives  some 

analytical  results  based  on  the  LSJ  process. 

2.  The  Gauss-Markov  Error  Process 

The  Gauss-Markov  process,  u(t),  introduces  explicit 

dependence  of  the  value  of  y(t+s)  based  on  the  value  of  y(t) 

for  any  s.  The  process  is  Markovian  in  that  the  present 

value  depends  only  on  the  value  that  immediately  proceeds 

it.  The  process  is  Gaussian  in  that,  for  any  n,  the  joint 

distribution  of  (f;(ti)},  i=l,...,n  is  multi variate-normal 

2 

each  with  zero  mean  and  variance  a  .  In  simulating  the 
process,  we  can  use  the  functional  form: 

y(t+s)  ■  e  XsU(t)  +  kn  (1.8) 

where :  1 

1.  k  -  (1  -  e“2Xs  )2  (see  Appendix  E) 

2 .  n  is  a  noripal  random  variable  with  zero  mean  and 
variance  a2. 
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The  process  results  in  a  continuous  sample  path,  but 
is,  in  fact,  nowhere  differentiable.  The  desirability  of 
using  this  type  of  error  function  in  modeling  acoustic  fluc¬ 
tuations  is  that  it  is  "more  like  reality".  Whether  it  is  in 
any  sense  "better"  than  the  LSJ  process  has  not  been 
determined.  Results  of  other  simulations  indicate  that 
cumulative  probabilities  of  detection  are  higher,  and  mean 
time  to  detect  is  lower,  when  the  GMA  process  is  used  instead 
of  the  LSJ.  However,  the  differences  are  not  significant  to 
the  point  of  determining  which  process  is  most  appropriate 
as  an  acoustic  model.  One  drawback  of  the  model  is  that  it 
is  not  possible  to  exactly  duplicate  a  continuous  function 
on  a  digital  computer,  and  thus  the  machine  represen tat ie- 
of  the  GMA  is  only  a  close  approximation  of  the  continuous 
error  process.  Thus,  the  manner  in  which  the  GMA  process  is 
digitally  replicated  has  an  effect  on  the  results  of  the 
simulation. 

The  parameterization  of  the  process  is  exactly 
analogous  to  that  of  the  LSJ  process.  Again,  X^  ^  can  be 
thought  of  as  rate  parameters,  and  a  ^  *  can  be  thought  of 
as  scale  parameters,  exactly  as  in  the  LSJ  model.  In  some 
analyses,  1/X  is  referred  to  as  the  "relaxation  time”  (and 
X  as  the  "relaxation  coefficient”)  since  it  measures  the 
decay  time  of  statistical  dependence  between  random  variables 
in  the  fluctuation  process  for  both  the  GMA  and  LSJ  processes 


The  covariance  function  for  the  Gauss-Markov  process 


is  exactly  the  sane  as  that  for  the  LSJ  (see  Equation  1.5) , 
and  therefore,  the  GMA  is  also  a  second  order  stationary 
process . 

There  are  no  analytical  results  similar  to  those 
from  the  LSJ  process  when  the  Gauss-Markov  process  is  used. 

See  Appendix  E  for  a  description  of  the  computer  representa¬ 
tion  of  the  GMA  process. 

G.  A  THREE-OUT-OF-FIVE  DETECTION  CRITERIA  MODEL 

McCabe  [Ref.  10:p.50]  contends  that  when  a  target  is 
exposed  for  only  a  short  amount  of  time  but  at  a  relatively 
high  signal-to-noise  ratio  (SNR)  ,  the  threshold  crossing  model 
returns  an  unrealistically  high  probability  of  detection,  and 
presents  a  signal  integration  model  based  on  a  time  dependent 
recognition  differential.  PASS  has,  as  an  option,  a  much 
simplified  model  which  is  easier,  and  faster  to  implement. 

It  seems  necessary  to  alter  the  threshold  crossing  model 
under  the  above  conditions  so  that  detection  would  require 
some  minimum  combination  of  signal  strength  and  time.  This, 
in  effect,  would  model  what  is  actually  integration  time  for 
the  sonar  processor.  In  general,  the  desired  relationship 
is  that  the  stronger  a  signal  is,  the  shorter  the  time 
required  for  recognition  as  a  valid  contact.  This  model 
seems  to  be  especially  necessary  when  detection  is  based  on 
other  them  aural  cues,  as  is  the  case  in  modern  sonars. 


As  a  user  option,  instead  of  the  threshold  crossing 
model  (where  detection  occurs  the  first  time  SE(t)  is  equal 
to  or  greater  than  zero) ,  a  3-of-5  model  can  be  employed 
which  requires  SE(t)  ^  0  on  3  of  the  last  5  samples  of 
signal  level  for  detection.  This  Minimum  Signal  Excess 
Logic  Model  (MSEL)  results  in  the  desired  effect  of  a  weak 
signal  being  required  to  be  "present"  for  a  longer  period 
for  detection  to  occur.  See  Appendix  F  for  a  more  complete 
analysis  of  the  model. 

McCabe  [Ref.  10:p.50]  recommends  that  an  integration 
model  be  used  whenever  strong  convergence  zones  are  present. 

H.  A  SEARCH  MEASURE  OF  EFFECTIVENESS  (MOE) 

1.  Searcher  MOE 

We  may  wish  to  measure  the  effectiveness  of  search 
(or  evasion)  in  terms  other  than  the  probability  of  secure 
detection.  For  example,  we  may  wish  to  penalize  the  searcher 
heavily  for  a  secure  target  counter-detection,  less  heavily 
for  a  simultaneous  detection,  and  still  less  heavily  for  no 
detection  (escape) . 

Let: 

1.  Ns  —  the  number  of  secure  searcher  detections 

2.  Nt  =  number  of  secure  target  detections 

3.  Nb  *  the  number  of  simultaneous  detections 

4.  Nn  *  number  of  no  detection  replications 

5.  Nr  *  the  total  number  of  simulation  replications 


Then  define  a  MQE  for  the  searcher  as: 


MOEs  = 


Ns-3Nt-2Nb-Nm  , 

- 3n? - + 


0.75 


where 


(1.9) 


OSMOEsSl  and  MOEs  -  1 Ns  *  Nr 

Note  that  the  weighting  factors  are  a  function  of  the 
relative  penalties  assigned  to  situations  other  than  a 
secure  detection,  and  could  be  easily  adjusted. 

2.  Target  MQE 

In  the  area  search  scenario,  we  assume  the  target  to 
be  aggressive,  and  therefore  a  target  MOE  is  assigned  analo¬ 
gous  to  the  searcher  MOE: 

M0Et  *  Nt~3N4ftf— ~^m  +  0.75  (1.10 

In  the  barrier  search,  the  target  is  assumed  to  want 
to  avoid  contact  with  the  searcher,  so  the  target  MOE  is 
adjusted  as  follows: 


™  _  Nt-3Ns-2Nb+Nm  ,  „  „ 

MOEt  =  - JUT -  +  °*75 

where,  in  both  cases. 


(1- 111 


N. 


-  Nr 


0<MOEtSl  and  MOEt 


1 


3.  The  Exchange  Ratio 


The  exchange  ratio  is  defined  as: 

MOE 

ER  *  (1-12) 

which  provides  a  composite  MOE  for  the  search. 

In  area  search: 

ER>1  —  Ns >Nt 
ER  =  1  <*>  Ns=Nt 
ER  Ns<Nt 

In  barrier  search: 

ER>1  ++■  Ns>Nt+jNn 
ER  *  1  NS»Nt+jNn 
ER<1  NS<Nt+jNn 

I.  PROGRAM  RESULTS /OUTPUT  AND  STATISTICAL  ANALYSIS 
The  following  results  are  available  from  PASS : 

1.  A  complete  record  of  the  input  data. 

2.  Probability  of  detection  by  the  searcher 

a.  Fraction  of  detections  which  were  direct  path. 

b.  Fraction  of  detections  which  were  in  the 
convergence  zones. 

c.  Fraction  of  detections  in  each  convergence  zone. 

3.  Probability  of  counter-detection  by  the  target. 

a.  Fraction  of  counter-detections  which  were  DP. 

b.  Fraction  of  counter-detections  which  were  in  the 
CZs. 


c.  Fraction  of  counter-detections  in  each  CZ. 

4.  Histogram  and  statistics  of  time  to  detection  by  searcher 

5.  Histogram  and  statistics  of  time  to  counter-detection 
by  the  target. 

6.  Histogram  and  statistics  of  range  at  detection  by 
searcher. 

7.  Histogram  and  statistics  of  range  at  counter-detection 
by  target. 

8.  Sectioning  of  time  and  range  data  (if  sample  size  is 
large  enough)  to  allow  for  estimates  of  the  variability 
of  sample  statistics. 

9.  A  plot  of  cumulative  probability  of  detection  versus  time 

10.  A  plot  of  cumulative  probability  of  counter-detection 
versus  time. 

11.  A  plot  of  conditional  probability  of  detection  versus 
range. 

12.  A  plot  of  the  searcher  positions  when  the  searcher 
detects  the  target. 

13.  A  plot  of  the  target  positions  when  the  searcher 
detects  the  target. 

14.  A  plot  of  conditional  probability  of  counter-detection 
versus  range. 

15.  A  plot  of  target  positions  when  target  counter-detects 
searcher. 

16.  A  plot  of  the  searcher  positions  when  the  target 
counter-detects  the  searcher. 

17.  A  plot  of  typical  searcher  and  target  FOM  for  a  five 
hour  period. 


J.  VALIDATION  OF  THE  PROGRAM 

There  is  no  existing  body  of  data  or  analysis  that  would 
lend  itself  to  a  direct  comparison  to  PASS  results  for  a 
rigorous  validity  check.  Numerous  area  and  barrier  search 


scenarios,  with  various  environmental,  platform,  and 
"random"  parameters,  were  run  with  results  not  inconsistent 
with  operational  experience,  limiting  theoretical  results, 
and  other  simulation  program  results.  Due  to  the  large 
number  of  variations  the  program  allows,  not  all  could  be 
tested.  The  strongest  statement  that  can  be  made  regarding 
PASS  validity  is  that  it  has  not  shown  to  be  invalid  after 
running  in  various  configurations. 


K.  PROGRAM  RUN  TIME 

PASS  is  a  relatively  expensive  analysis  tool  from  the 
standpoint  of  computer  run  time.  On  the  NPGS  computer  (IBM 
VM-370  with  three  3033  CPU's)  a  rule  of  thumb  is  that  it 
will  take  about  50  CPU-seconds  per  1000  replications  of  an 
area  search  in  100  nm  x  100  nm  search  area.  Obviously,  more 
or  less  time  will  be  required  for  different  sensors  and 
other  scenarios.  The  point  to  be  made  is  that  the  run  time 
is  long  enough  to  warrant  extensive  pre-planning  so  that  the 
maximum  desired  information  is  obtained  from  each  run. 


II.  DESCRIPTION  OF  PROGRAM  PASS 

A.  INTRODUCTION 

The  basic  flow  diagram  for  program  PASS  is  shown  in 
Figures  2.1,  2.2,  and  2.3.  Each  "part"  of  the  program,  as 
described  in  the  comment  blocks  in  the  source  code  (see 
Appendix  H)  is  indicated  in  parentheses  in  each  block  of  the 
flow  diagram.  The  flow  diagram,  the  annotated  source  code, 
and  the  description  of  the  variable  names  (see  Appendix  B) 
should  be  used  with  this  chapter  to  gain  an  understanding  of 
the  logic  employed  in  the  simulation  and  the  implementation 
of  this  losic  in  the  FORTRAN  code. 

B.  DATA  I$PUT  AND  SIMULATION  INITIALIZATION  (PROGRAM 

PARTS  1,  2,  3) 

All  logical  variables  are  set  to  default  values  prior  to 
data  input. 

Data  input  is  normally  accomplished  by  the  user  respond¬ 
ing  to  a  number  of  options  presented  on  the  terminal  screen. 
These  options  are  displayed  by  sequential  calls  to  sub¬ 
routines  OPTNO  through  0PTN2G.  During  this  interactive 
session,  the  data  necessary  to  run  the  program  in  the 
default  mode  is  written  to  file  07  on  the  A-disk,  which 
allows  saving  the  basic  data  from  one  run  to  the  next.  The 
user  can  bypass  all  the  options  on  the  screen  and  proceed 
directly  to  program  execution  by  reading  the  data  from  this 
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file  through  the  subroutine  REAOIT.  Invoking  subroutine 
READIT  requires  accepting  the  following  defaults: 

1.  A  complete,  and  properly  formatted,  input  data  set  is 
assumed  to  be  on  the  A-disk  in  file  PASS  DATA  (file 
07)  . 

2.  All  input  data  will  be  sent  to  the  output  file  PASS 
OUTPUT  (file  06) . 

3.  The  acoustic  fluctuation  model  is  Lambda-Sigma  Jump 
(LSJ) . 

4.  A  complete  statistical  analysis  of  the  results  is  to 
be  done. 

5.  The  3-of-5  detection  criteria  model  is  not  used. 

6.  For  area  search  the  initial  target  position  is 
distributed  uniform  on  (0,XMAX) ,  (0,YMAX). 

7.  For  barrier  search,  the  initial  target  lateral  position 
is  distributed  uniform  on  (Q,XMAX). 

8.  The  searcher  figure-of -merit  is  independent  of  target 
aspect . 

Subroutines  ECHOl  through  ECH06  write  the  input  data  to 
the  output  file  if  the  user  so  desires. 

The  limits  of  the  target  course  variations  in  the 
barrier  scenario  (ANG1,  ANG2)  are  initialized  if  barrier 
search  is  selected.  Detection  counters  (NDO,  NDT,  NCZDO(I) , 
NCZDT(I) ,  NBOTH,  NONE)  are  set  to  zero,  and  the  stack  pointers 
(Cl,  C2,  C3,  C4)  are  initialized.  The  direction  of  each 
searcher  path  leg  (DX(I)  ,  DY(I))  and  the  distance  of  each 
leg  (DIST(I) )  are  computed,  and  the  target  speed  range 
(STINC)  and  boundary  reflection  constants  (TOXMAX,  TOYMAX) 
are  calculated. 
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C.  REPLICATION  COUNTER  AND  SETUP  (PROGRAM  PART  4) 

If  the  required  number  of  replications  (NREP)  have  been 
completed,  transfer  is  made  to  the  data  analysis  and  output 
routines  (Part  18)  and  subsequent  program  termination. 
Otherwise,  another  replication  is  started. 

A  count  of  completed  replications  is  written  to  the 
screen  every  200  replications.  Counters  for  the  collection 
of  a  representative  sample  of  figure-of -merit  data„„iMi";""M2 ) 
are  reset  if  five  hours  of  this  data  harsfnot  been  collected 
in  a  previous  replication.  Simulation  times  (TNOW,  TLAST) 
are  set  to  zero,  and  signal  MSEL  counters  (MSELO,  MSELT)  and 
convergence  zone  (CZ)  detection  pointers  (KCZO,  KCZT)  are 
reset . 

D.  INITIALIZATION  OF  THE  REPLICATION  (PROGRAM  PART  5) 

At  the  start  of  each  replication,  the  following  initial¬ 
ization  and  setup  is  accomplished: 

1.  The  time  to  the  next  fluctuation  change  (TIFL(I) , 
1-1,2, 3)  is  selected  from  an  exponential  distribution, 
parameterized  by  ALAM(I) ,  by  calls  to  sub routine  EXPO. 

2.  The  initial  magnitude  of  the  fluctuation  (AFL(I))  is 
selected  from  a  normal  distribution,  parameterized  by 
SIGMA(I) ,  by  calls  to  the  subroutine  XLS. 

3.  The  searcher  position  (XO,  YO)  is  set  to  the  first 
searcher  track  anchor  point  (XP(1),  YP(1)). 

4.  The  search  leg  pointer  (NLEG)  is  set  to  one,  and  the 
search  mode  (MODE)  is  set  to'  "drift”.  The  time  to  the 
next  searcher  speed  change  (TSC)  is  set  to  the  drift 
time.  The  figure-of-merit  for  both  platforms  (FOMO, 
FOMT)  is  selected  based  on  searcher  drift  speed,  and 
searcher  speed  (SO)  is  set  to  drift  speed. 


5.  The  tine  to  the  searcher  course  change  (TCC)  is  set 
based  on  the  first  leg  distance  and  drift  speed,  and 
searcher  speed  vectors  (VX,  VY)  are  set  based  on  the 
direction  vectors  and  search  speed. 

6.  If  the  user  is  to  specify  the  initial  target  lateral 
position  in  the  barrier  scenario,  the  initial  target 
X-position  (XT)  is  selected  by  a  call  to  subroutine 
XDISTB.  Otherwise,  the  initial  target  X-position  is 
selected  uniform  on  (0,XMAX)  by  a  call  to  subroutine 
UZ1. 

7.  If  the  barrier  scenario  is  selected,  the  initial 
target  Y-position  (YT)  is  selected  a  distance  offset 
from  YMAX.  This  offset  distance  is  representative  of 
the  searcher  being  "tine-late"  on  the  barrier.  The 
offset  is  calculated  based  on  the  maximum  time  late 
(START) ,  the  mean  target  speed,  and  a  call  to  sub¬ 
routine  UZ1.  In  the  area  search  scenario,  the  initial 
target  Y-position  is  selected  uniform  on  (0,YMAX)  by 

a  call  to  subroutine  UZ1. 

8.  Target  speed  (ST)  is  selected  uniform  on  (STMIN,STMAX) 
by  a  call  to  subroutine  UZ1. 

9.  The  time  to  target  speed  change  (TTSC)  is  selected 
from  the  exponential  distribution,  parameterized  by 
RTSC,  by  a  call  to  subroutine  EXPO. 

10.  In  the  barrier  scenario,  the  initial  target  psuedo- 
course  (X)  is  selected  uniform  on  (ANGl ,ANG2) ,  and  in 
the  area  search  scenario,  uniform  on  (0,360),  by  a 
call  to  subroutine  UZ1.  If  the  searcher  FOM  is 
dependent  on  target  aspect ,  the  target  psuedo-course  is 
saved  (THETA) .  Target  direction  vectors  (COSX,  SINX) 
and  speed  vectors  (UX,  UY)  are  computed  based  on 
target  course  and  speed. 

11.  Uxe  time  to  target  course  change  (TTCC)  is  selected 
from  an  exponential  distribution,  parameterized  by 
RTCC,  by  a  call  to  subroutine  EXPO. 


E.  DETERMINE  DETECTION  RANGES  (PROGRAM  PART  6) 

If  the  searcher  FOM  is  dependent  on  target  aspect,  the 
searcher  FOM  computed  in  Part  5  or  14  is  re-calculated  based 
on  the  current  relative  bearing  from  the  target  to  the 
searcher  (BREL) .  The  subroutine  RELB  returs  BREL  based  on 


36 


the  current  searcher-target  geometry.  The  subroutine  INTRPL 
returns  interpolated  values  of  searcher  FQM  (FOMO)  based  on 
BREL  and  the  searcher-speed/target-relative-bearing  FOM  data 
(FOMBD(I) ,  FOMBS(I)).  The  searcher  FOM  (X)  is  computed  by 
applying  the  error  process  values  (AFL(l) ,  AFL ( 2 ) )  to  FOMO. 
Five  hours  of  searcher  FOM  (FO(I))  is  saved  for  graphical 
display.  The  direct-path  (DP)  detection  range  for  the 
searcher  (RNGO)  is  obtained  by  applying  the  searcher  FOM  (x) 
to  the  DP  propagation  loss  curve  through  subroutine  INTRPL. 

If  convergence  zones  are  present  (NCZO  not  zero) ,  the 
searcher  FOM  (X)  is  applied  to  each  CZ  propagation  loss  value 
(CZLO(I))  to  determine  the  most  distant  CZ  in  which  a  detec¬ 
tion  can  take  place  (KCZO) .  If  KCZO  is  zero,  then  no 
searcher  detection  can  take  place  in  a  CZ. 

The  preceeding  procedure,  with  the  exception  of-  the 
aspect  dependency  portion,  is  applied  also  to  the  target 
related  data  to  produce  values  for  RNGT  and  KCZT. 

Hie  maximum  range  at  which  a  direct-path  detection  can 
take  place  by  either  platform  (RMAX)  is  calculated  from  the 
maximum  of  RNGO,  RNGT. 

F.  CHECK  FOR  DETECTION  (PROGRAM  PART  7) 

The  range  between  the  platforms  (RNG)  is  calculated 
based  on  the  current  geometry,  and  the  algebraic  difference 
between  this  range  and  the  maximum  DP  detection  range  is 
calculated  (TBIG) .  Direct-path  detection  is  possible  for  at 
least  one  platform  only  if  TBIG  is  less  than  or  equal  to  zero 


If  the  searcher  FQM  will  support  detection  in  one  or 
more  CZs  (KCZO  not  zero) ,  the  range  (RNG)  is  checked  to  see 
if  the  target  is  in  any  searcher  CZ  annulus,  starting  from 
the  most  distant  one  that  will  support  detection,  and 
working  in  toward  the  searcher.  If  the  target  is  in  one  of 
the  CZs,  KCZO  is  set  to  minus-one  (as  a  CZ  detection  flag) , 
and  the  CZ  in  which  detection  is  possible  (ICZO)  is  saved. 

An  identical  CZ  detection  check  for  the  target  is  done, 
and  if  the  target  can  detect  the  searcher  in  a  CZ,  KCZT  is 
set  to  minus  one,  and  ICZT  is  saved. 

The  value  of  T3IG  is  set  to  the  minimum  distance  the  plat¬ 
forms  would  have  to  traverse  (assuming  a  head-to-head  closure) 
for  a  detection  to  just  occur,  based  on  the  current  searcher 
and  target  speeds  and  FOMs.  That  is,  if  the  searcher  and 
target  were  headed  directly  toward  each  other,  TBIG  would  be 
the  total  distance  covered  by  both  platforms  at  the  time  of 
first  detection  by  either. 

If  a  CZ  detection  is  possible  (KCZO  and/or  KCZT  less 
than  zero)  or  if  a  DP  detection  is  possible  (TBIG  less  than 
zero) ,  a  transfer  is  made  to  Part  16A  where  it  is  determined 
if  a  detection  is  made,  and  by  whom. 

If  detection  is  not  possible,  transfer  is  made  to 
Part  8. 

G.  COMPOTE  THE  TIME  OF  THE  "NEXT  EVENT"  (PROGRAM  PART  8) 

If  the  MSEL  model  is  used,  and  either  MSEL  counter  (MSELO, 
MSELT)  is  not  in  a  zero  state,  this  implies  that  detection 
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was  possible  at  some  earlier  time,  and  now  it  is  not. 
Therefore,  any  non-zero  MSEL  counter  is  decremented  by  one. 
If,  after  decrementing,  any  MSEL  counter  is  not  in  a  zero 
state,  TBIG  is  set  to  zero.  This  results  in  the  maximum 
simulation  time-step  defaulting  to  0.05  hours. 

The  maximum  simulation  time-step  is  calculated  by  taking 
the  maximum  of  0.05  hours  (3  minutes)  or  the  time  required 
for  searcher  and  target  to  close  the  distance  TBIG  assuming 
all  of  the  current  platform  speeds  (SO,  ST)  are  in  the  line 
of  sight  (head-to-head  closure) .  The  range  TBIG  is  replaced 
by  the  time  TBIG.  The  incremented  time  (TINC)  is  set  to  the 
current  time  (TNOW)  plus  TBIG,  and  the  next  event  pointer 
(J)  is  set  to  8.  Next,  the  times  of  all  other  possible 
events  are  checked  to  see  if  any  are  earlier  than  TINC.  If 
an  event  time  is  earlier  than  TINC,  the  time  of  this  event 
replaces  the  current  value  of  TINC,  and  the  next  event 
pointer  is  changed  appropriately.  When  this  process  is 
completed,  TINC  will  be  the  time  of  the  "next  event",  and 
the  value  of  J  will  identify  what  that  event  is.  Note  that 
if  J  remains  at  8,  the  time  increment  will  be  TBIG.  TBIG 
tends  to  be  large  when  the  platforms  are  far  apart,  and 
decreases  as  the  range  decreases. 

If  TINC  exceeds  the  maximum  allowable  search  time  (TMAX) , 
then  transfer  is  made  to  Part  17  for  the  recording  of  no¬ 
detection  data,  and  the  start  of  a  new  replication,  if 
appropriate . 


H.  MOVE  THE  SHIPS  (PROGRAM  PARTS  9A,  9B) 

The  time  increment  (X)  is  obtained  by  subtracting  the 
current  time  from  the  incremented  time,  and  the  projected 
searcher  position  (XOT,  YOT)  is  obtained  by  applying  the 
time  increment  to  the  searcher  velocity  vectors  (VX,  VY) . 

The  distance  from  the  current  searcher  position  (XO,  YO)  to 
the  projected  searcher  position  (Dl)  is  calculated  from  the 
geometry.  The  distance  from  the  current  searcher  position 
to  the  next  searcher  track  anchor  point  (D2)  is  calculated 
from  the  geometry.  If  Dl  is  less  than  or  equal  to  D2,  the 
searcher  position  is  updated  to  the  projected  position.  If 
Dl  is  greater  than  D2,  then  the  time  increment  (X)  is  adjusted 
to  put  the  searcher  at  the  next  anchor  point,  and  the  next 
event  pointer  (J)  is  set  to  4  which  will,  in  Part  12,  cause 
the  searcher  position  to  be  updated  to  the  anchor  point. 

The  target  position  (XT,  YT)  is  updated  by  applying  the 
time  increment  (X)  to  the  target  velocity  vectors  (UX,  UY) . 

The  simulation  time  (TNOW)  is  updated  by  the  time 
increment  (X) . 

If  the  target  position,  as  calculated  in  Part  9A,  ends 
up  outside  the  search  area  (0,XMAX) ,  (0,YMAX),  the  target 
position  and  velocity  vectors  are  adjusted  such  that  the 
boundaries  appear  to  "reflect"  the  target.  The  target  track 
behaves  as  would  a  light  beam  striking  the  mirrored  surfaces 
of  the  search  area  boundaries.  If  the  searcher  FOM  is  target 


aspect  dependent,  then  the  new  target  course  must  be  saved 
(THETA)  for  use  in  subsequent  calculation  of  relative 
bearing. 


I.  BRANCH  TO  THE  "NEXT  EVENT"  (PROGRAM  PART  10) 

Depending  on  the  value  of  the  next  event  pointer  ( J) , 
trams fer  is  made  to  the  program  part  which  accomplishes  the 
physical  event  associated  with  the  time-step  computed  in 
Part  8. 

If  J»8  (which  means  "no  event”  takes  place) ,  transfer  is 
made  to  Part  6  or  Part  7,  depending  on  the  integration  model 
in  effect.  If  the  MSEL  (3-of-5)  model  is  in  effect, 
transfer  must  be  made  to  Part  6  to  reset  the  KCZO  and  KCZT 
flags.  In  future  sections  where  transfer  is  made  to  Part  6 
or  7,  the  same  dependence  on  the  MSEL  detection  model  in 
effect  will  govern  where  transfer  is  made. 

If  J»l,2,3  then  2m  acoustic  fluctuation  level  is  to 
change,  and  transfer  is  made  to  Part  11. 

If  J»4  then  the  searcher  course  is  to  change  at  a 
searcher  track  anchor  point,  and  transfer  is  made  to  Part  12 

If  J»5  then  the  target  course  is  to  change,  and  transfer 
is  made  to  Part  13. 

If  J»6  then  the  target  speed  is  to  change,  and  transfer 
is  made  to  Part  15. 

If  J«7  then  the  searcher  speed  is  to  change,  and 
transfer  is  made  to  Part  14. 
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J.  ACOUSTIC  FLUCTUATION  LEVEL  CHANGE  (PROGRAM  PART  11) 

The  fluctuation  level  change  is  calculated  differently 

depending  on  the  model  used.  The  next  event  pointer  (J)  also 
identifies  which  error  term  is  to  change. 

In  the  LSJ  model,  a  call  to  the  subroutine  XLS  with 
input  parameter  SIGMA (J)  returns  a  new  fluctuation  level 
(AFL(J) )  from  a  normal  distribution.  Next,  a  call  to  the 
subroutine  EXPO  with  input  parameter  ALAM(J)  returns  a  new 
time  increment  (T)  to  the  next  fluctuation,  which  is  added 
to  current  time  to  obtain  the  time  ofthe  next  fluctuation 
(TIFI(J)).  Transfer  is  then  made  to  Part  6. 

In  the  GMA  model,  the  time  to  the  next  fluctuation  is 
obtained  in  a  manner  identical  to  the  LSJ  model.  The  lag 
time  from  the  last  calculation  of  error  signals  (S)  is 
obtained  by  subtracting  the  last  time  the  levels  were  calcu¬ 
lated  (TLAST)  from  the  current  time,  and  then  setting  TLAST 
to  the  current  time.  All  three  of  the  signal  fluctuations 
are  calculated  by  successive  calls  to  subroutine  XLS  parame¬ 
terized  by  SIGMA(JJ) ,  JJ=1,2,3,  and  the  algebraic  expression 
developed  in  Appendix  E  which  involves  ALAM(JJ) ,  JJ=1,2,3. 
Trams fer  is  then  made  to  Part  6. 

K.  SEARCHER  COURSE  CHANGE  (PROGRAM  PART  12) 

The  searcher  position  (XO,  YO)  is  assigned  to  the  next 
search  leg  start  point  (XP(NLEG) ,  YP(NLEG))  and  new  searcher 
velocity  vectors  (VX,  VY)  are  calculated  based  on  the  next 
leg  direction  vectors  (DX(NLEG),  DY(NLEG))  and  current 


search  speed  (SO) .  The  time  to  the  next  searcher  course 
change  (TCC)  is  calculated  based  on  the  leg  distance 
(DIST(NLEG) )  and  current  searcher  speed.  Transfer  is  made 
to  Part  6  ot  7,  depending  on  the  detection  model  in  effect. 

L.  TARGET  COURSE  CHANGE  (PROGRAM  PART  13) 

If  the  area  search  scenario  is  in  effect,  the  new  target 
course  (X)  is  selected  from  a  uniform  distribution  on 
(0,360)  by  a  call  to  subroutine  UZ1.  If  the  barrier  scenario 
is  in  effect,  the  new  target  course  is  selected  from  a 
uniform  distribution  on  (ANG1,ANG2) .  The  new  target  course 
is  saved  (THETA)  if  the  searcher  FOM  is  dependent  on  target 
aspect . 

New  target  speed  vectors  (UX,  UY)  are  calculated  based 
on  the  new  course  and  current  target  speed  (ST) .  The  time 
interval  to  the  next  target  course  change  (T)  is  obtained 
from  an  exponential  distribution  by  calling  the  subroutine 
EXPO  with  input  parameter  RTCC,  and  the  time  of  the  next 
course  change  (TTCC)  is  calculated  based  on  this  interval. 

If  the  searcher  FOM  is  a  function  of  target  aspect, 
transfer  is  made  to  Part  6  (via  Part  11  to  calculate  new 
fluctuation  levels  if  GMA  is  the  fluctuation  model)  to 
determine  new  detection  ranges.  Otherwise,  transfer  is 
made  to  Part  6  or  7  depending  on  the  detection  model  in 


effect 


M.  SEARCHER  SPEED  CHANGE  (PROGRAM  PART  14} 

The  MODE  state  is  changed  to  reflect  the  new  searcher 
speed.  The  searcher  and  target  FOM  ( FOMO ,  FOMT)  are  changed 
to  reflect  the  new  searcher  speed,  and  the  search  speed  (SO) 
is  updated.  The  tine  of  the  next  searcher  speed  change  (TSC) 
is  updated  using  the  equivalence  variable  TIME (MODE) .  The 
time  to  the  next  searcher  course  change  is  updated  based  on 
the  new  search  speed,  and  distance  to  the  next  searcher 
track  anchor  point  (D) .  Transfer  is  made  to  Part  6  (via 
Part  11  if  GMA  is  the  fluctuation  model)  to  determine  new 
detection  ranges. 

N.  TARGET  SPEED  CHANGE  (PROGRAM  PART  15) 

The  new  target  speed  (ST)  is  selected  from  a  uniform 
distribution  on  (STMIN,  STMAX)  by  a  call  to  the  subroutine 
UZ1.  The  time  interval  to  the  next  speed  change  (T)  is 
selected  from  an  exponential  distribution  by  a  call  to  the 
subroutine  EXPO  with  input  parameter  RTSC,  and  the  time  of 
the  next  speed  change  (TTSC)  is  calculated  based  on  this 
interval.  New  target  speed  vectors  (UX,  UY)  are  calculated 
based  on  the  new  speed  and  current  target  direction  vectors 
(COSX,  SINX) .  Transfer  is  then  made  to  Part  6  or  7, 
depending  on  the  detection  model  in  effect. 

O.  DETERMINE  WHICH  PLATFORM  CAN  DETECT  (PROGRAM  PART  16A) 

This  part  of  the  program  is  entered  only  if,  in  Part  7, 

detection  by  at  least  one  of  the  platforms  is  possible. 


Therefore,  if  the  direct  path  detection  range  of  the  target 
(RNGT)  is  less  than  the  current  range  (RNG) ,  and  the  target 
cannot  make  a  CZ  detection  on  the  searcher  (KCZT  greater  or 
equal  to  zero) ,  the  detection  opportunity  must  exclusively 
belong  to  the  searcher,  and  transfer  is  made  to  Part  16C. 
Similarly,  if  the  direct  path  detection  range  of  the 
searcher  CRNGO)  is  less  them  the  current  range,  and  the 
searcher  cannot  make  a  CZ  detection  on  the  target  (KCZO 
greater  or  equal  to  zero) ,  the  detection  opportunity  must 
exclusively  belong  to  the  target,  and  transfer  is  made  to 
Part  16D. 

If  neither  of  the  above  cases  are  true,  then  the  detec¬ 
tion  must  be  simultaneous,  and  trans fer  is  made  to  Part  16B. 

P.  SIMULTANEOUS  DETECTION  OPPORTUNITY  (PROGRAM  PART  16B) 

If  MSEL  is  not  in  effect,  the  simultaneous  detection 
counter  (NBOTH)  is  incremented,  and  transfer  is  made  to 
Part  4  to  start  a  new  replication  or  output  results. 

If  the  MSEL  model  is  used,  both  the  searcher  and 
target  counters  (MSELO,  MSELT)  are  incremented.  If  both 
the  counters  equal  three,  then  the  actions  in  the  preceeding 
paragraph  are  carried  out.  If  only  the  searcher  counter  is 
at  three,  then  transfer  is  made  to  Part  16C.  If  only  the 
target  counter  is  at  three,  transfer  is  made  to  Part  16D. 

If  no  -  detection  takes  place  (both  counters  less  than  3) 
transfer  is  made  to  Part  8  to  determine  the  time-step  to 


the  next  event 


SEARCHER  SECURE  DETECTION  OPPORTUNITY  (.PROGRAM  PART  16C) 
If  the  MSEL  model  is  not  used,  the  following  takes  place 

1.  Searcher  detection  counter  (NDO)  is  incremented. 

2.  If  this  was  a  CZ  detection  (KCZO  *  -1) ,  the  number  of 
detections  in  the  particular  CZ  is  incremented 
(NCZDO(ICZO) ) . 

3.  The  time  of  the  detection  is  saved  (NTDO (NDO) ) . 

4.  The  range  of  the  detection  is  saved  (RNTDO(NDO) ) . 

5.  Searcher  position  at  time  of  detection  is  saved 
(XODT(NDO) ,  YODT(NDO) ) . 

6.  Target  position  at  time  of  detection  is  saved 
(XTOD (NDO) ,  YTOD (NDO) ) . 

7.  Transfer  is  then  made  to  Part  4  to  begin  a  new 
replication  or  output  results. 

If  the  MSEL  model  is  used,  the  following  takes  place: 

1.  Searcher  MSEL  counter  (MSELO)  is  incremented. 

2.  If  MSELO  is  3,  the  actions  in  the  preceeding  paragraph 
are  executed. 

3.  If  MSELO  is  less  than  three,  the  target  MSEL  counter 
(MSELT)  is  decremented  if  it  was  not  zero,  and 
transfer  is  made  to  Part  8. 


TARGET  SECURE  DETECTION  OPPORTUNITY  (PROGRAM  PART  16D) 

If  the  MSEL  model  is  not  used,  the  following  takes  place 

1.  Target  detection  counter  (NDT)  is  incremented. 

2.  If  this  was  a  CZ  detection  (KCZT  =  -1) ,  the  number  of 
detections  in  the  particular  CZ  is  incremented 
(NCZDT(ICZT) ) . 

3.  The  time  of  the  detection  is  saved  (NTDT(NDT)). 

4.  The  range  of  the  detection  is  saved  (RNTDT(NDT) ) . 

5.  Searcher  position  at  time  of  detection  is  saved 
(XTDO(NDT) ,  YTDO(NDT) ) . 


6.  Target  position  at  time  of  detection  is  saved 
(XOTD(NDT) ,  YOTD(NDT) ) . 

7.  Transfer  is  then  made  to  Part  4  to  begin  a  new 
replication  or  output  results. 

If  the  MSEL  model  is  used,  the  following  takes  place: 

1.  Target  MSEL  counter  (MSELT)  is  incremented. 

2.  If  MSELT  is  3,  the  actions  in  the  preceeding  para¬ 
graph  are  executed. 

3.  If  MSELT  is  less  than  three,  the  searcher  MSEL 
counter  (MSELO)  is  decremented  if  it  was  not  zero, 
and  transfer  is  made  to  Part  8. 


S.  NO  DETECTION  OCCURS  (PROGRAM  PART  17) 

If  no  detection  occurs  by  the  maximum  allowed  search 
time  (TMAX) ,  the  no-detection  counter  (NONE)  is  incremented, 
and  transfer  is  made  to  Part  4  to  start  a  new  replication 
or  output  results. 


T.  RESULTS  OUTPUT  (PROGRAM  PART  18) 

By  calling  subroutine  SINKEM,  the  following  output 
functions  are  accomplished: 

1.  Some  results  of  the  simulation  (tiroes,  ranges)  are 
written  to  file  06  on  the  A-disk  in  the  form  of 
histograms  and  statistics . 

2.  Some  results  of  the  simulation  (times,  ranges, 
cumulative  probability  of  detection/counter¬ 
detection,  geographical  positions,  and  figures-of- 
merit)  are  written  to  file  10  on  the  A-disk.  File  10 
is  designed  to  be  a  data  file  for  the  FORTRAN  program 
PASPLT  which  produces  graphical  output  by  invoking 
the  DISSPLA  graphics  system. 

3.  An  abbreviated  summary  of  the  simulation  results  is 
sent  to  the  terminal  screen  (file  08)  and  to  file  06 
on  the  A-disk. 


Figure  2.1.  Flow  Diagram  for  Program  PASS  Parts  1-7 


III.  EXAMPLE  SEARCH  PROBLEMS 


A.  INTRODUCTION 

In  this  chapter  a  set  of  search  scenarios  will  be  inves¬ 
tigated  to  demonstrate  the  operation  and  capabilities  of 
PASS.  In  each  scenario  (with  the  exception  of  the  definite 
range  law  examples) ,  environmental  and  performance  parameters 
are  representative  of  "real-world"  conditions,  but  they  do 
not  represent  any  specific  set  of  platforms  or  environmental 
conditions.  This  was  done  to  preclude  questions  regarding 
security  classification. 

The  purpose  of  this  chapter  is  demonstration.  As  such, 
the  problems  investigated  are  treated  in  what  may  be  called 
a  superficial  manner.  That  is,  the  questions  posed  are  not 
rigorously  analyzed,  and  conclusions  drawn  can  only  be 
applied  within  the  narrow  constructs  of  the  specific  scenario. 

With  one  exception,  the  output  from  PASS  (and  PASPLT)  as 
shown  in  the  figures  in  this  chapter  is  unaltered.  The 
exception  is  the  histogram  and  statistics  output  produced  by 
the  NONIMSL  subroutine  HISTGP,  which  is  produced  on  the 
computer  center  line  printer.  The  format  of  the  histograms 
and  statistical  output  is  oversize,  and  is  not  of  high 
enough  quality  to  be  photo-reduced  for  legible  inclusion  in 
this  thesis.  Selected  statistical  results  from  HISTGP  will 
be  provided  in  tabular  form  when  required. 


B.  AREA  SEARCH:  NOISE  QUIETING  VS.  IMPROVED  SENSOR 

PERFORMANCE 

1.  Scenario  Description 

A  target  with  a  relatively  large  acoustic  disadvan¬ 
tage  is  the  subject  of  systematic  area  search.  The  question 
examined  is  which  of  the  following  possibilities  would  be 
more  advantageous  to  the  target: 

1.  Decreasing  target  radiated  noise  by  3  db,  or 

2.  Increasing  the  FOM  of  the  target  sonar  by  3  db. 

2.  PASS  Input  Data 

Table  I  shows  the  input  data  for  the  base  case  (case 
A-I) ,  in  the  format  presented  in  the  PASS  output  file. 

3.  Environment  and  Search  Plan 

Figures  3.1  and  3.2  show  the  propagation  loss  models 

used  in  the  example.  Note  that  even  before  figure-of-merit 
is  calculated,  the  searcher  has  an  advantageous  situation  in 
that  his  propagation  loss  curve  predicts  longer  detection 
ranges  than  for  the  target  (for  the  same  FOM) .  This  occurs 
because  the  geometric  mean  frequency  of  the  searcher  sonar 
is  less  than  that  of  the  target.  Due  to  the  presense  of 
convergence  zones,  the  MSEL  model  is  used. 

Figure  3.3  shows  the  searcher  track,  with  track 
spacing  of  75  miles  (about  equal  to  the  range  of  the  second 
CZ) .  Hie  search  area  is  300  nm  square. 

4.  Platform  Parameters 

The  search  speed  is  established  as  constant  at  10 


knots  by  setting  sprint  and  drift  speeds  equal  to  10  knots, 


and  searcher  sprint  and  drift  FOM  equal  at  87  db.  Whenever 
constant  speed  search  is  to  be  simulated ,  the  times  at  sprint 
and  drift  speeds  should  be  set  to  some  large  number  (500  hours 
in  this  case)  to  preclude  unnecessary  program  steps. 

The  target  speed  is  uniformly  distributed  between  5 
and  10  knots,  with  changes  in  course  and  speed  taking  place 
on  the  average  of  twice  per  hour.  The  target  FOM  is  84  db  in 
the  base  case. 

5.  Simulation  Run  Parameters 

Two  thousand  replications  were  run,  and  a  maximum 
search  time  was  set  at  720  hours  (30  days) .  The  random 
number  seed  can  be  any  integer,  and  for  these  runs  was 
selected  as  the  time  of  program  execution. 

6.  Results  and  Conclusions 

Numerical  results  of  the  base  case  (Case  A- I)  and 
modified  scenarios  (Cases  A-II,  A- III,  and  A- IV)  are  shown 
in  Table  II.  The  modified  scenarios  reflect  the  following 


changes  to  the  basic  case: 

1.  Case  A-II:  The  target  platform  gets  a  3  db  gain  in 
sensor  performance,  which  results  in  the  target  and 
searcher  having  equal  FOM  at  87  db. 

2.  Case  A-III:  The  target  platform  reduces  its  radiated 
noise  by  3  db  which  results  in  the  target  and  searcher 
having  equal  FOM  at  84  db. 

;  Case  A- IV:  For  comparison  purposes,  the  target  is 

assumed  to  achieve  both  the  improvements  in  the  above 
cases,  resulting  in  a  target  FOM  of  87  db,  and  a 
searcher  FOM  of  84  db. 


The  simulation  results  indicate  that  Case  A-III  is 

•  i 

slightly  more  favorable  to  the  target  than  Case  A- 1 1  in  the 
following  respects: 

1.  The  cumulative  probability  of  (secure)  detection  by 
the  searcher  is  smaller. 

2.  A  larger  proportion  of  the  searcher  detections  are 
direct  path,  resulting  in  decreased  mean  detection 
range,  and  increased  probability  of  counter-detection. 

3.  Less  detections  take  place  in  the  most  distant  CZ. 

4.  The  change  in  the  exchange  ratio  in  A-III  (92.2)  is 
greater  than  that  for  A-II  (83.2). 

5.  The  mean  time  to  detection  is  longer. 

6.  The  mean  and  median  detection  ranges  are  shorter. 

The  conclusion  is  that  a  greater  expected  increase 
in  performance  will  be  realized  if  the  3  db  gain  is  obtained 
by  reducing  radiated  noise.  Because  of  the  better  perfor¬ 
mance  of  the  searcher  sonar,  the  average  slope  of  the 
propagation  loss  curve  (change  in  prop,  loss/change  in 
detection  range)  will  be  less  them  that  for  the  target  (i.e. 
a  "flatter"  curve).  Thus,  the  change  in  detection  range 
per  3  db  change  would  be  greater  for  the  searcher  than  the 
target.  Therefore,  reducing  the  searcher  effectiveness  (by 
quieting  the  target)  has  a  larger  effect  on  detection 
performance  them  increasing  the  target  effectiveness  (by 
improving  the  sensor)  has  on  counter-detection  perf ormance . 

For  comparison  purposes,  a  complete  set  of  graphical 
output  from  PASS  for  Case  A-I  is  shown  in  Figure  3.4  through 
3.8,  and  for  Case  A-III  in  Figures  3.9  through  3.17. 
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TABLE  II 


NUMERICAL  RESULTS  FOR  CASE  A  EXAMPLES 


Case  A- I 

Case  A-ZI 

Case  A— III 

Case  A- IV 

PD 

0.9725 

0.9375 

0.9210 

0.8490 

PDDP 

0.1496 

0.1515 

0.1900 

0.1914 

PDCZ 

0.8504 

0.8485 

0.8100 

0.8086 

PDCZ1 

0.1172 

0.0960 

0.1482 

0.0960 

PDCZ2 

0.7322 

0.7525 

0.6618 

0.7126 

HOES 

0.9790 

0.9502 

0.9370 

0.8732 

PCD 

0.0015 

0.0115 

0.0150 

0.0540 

PCDDP 

0.3333 

0.0000 

0.0000 

0.0093 

PCDCZ 

0.6667 

1.0000 

1.0000 

0.9907 

PCDCZ1 

0.6667 

1.0000 

1.0000 

0.9907 

PCDCZ2 

0.0000 

0.0000 

0.0000 

0.0000 

MCET 

0.0080 

0.0243 

0.0310 

0.0783 

ER 

122.4 

39.19 

30.20 

11.16 

(T)S 

33.7 

35.2 

41.8 

42.9 

(T)  t 

26.4 

17.8 

32.3 

33.9 

(R)s 

68.9 

69.7 

64.6 

67.7 

(R-50)  s 

75.6 

75.6 

73.1 

74.1 

(R)t 

33.0 

38.0 

37.9 

38.3 

(R-50)  t 

4.75 

39.1 

37.9 

39.2 

Notes :  1.  Dus  to  the  low  number  of  counter-detections  in  Cases  A- 1 

and  A- I I  13  and  23  respectively)  the  variance  of  the  target 
parameter  estimates  will  be  large.  Therefore,  comparison 
should  be  based  on  searcher  parameter  estimates. 


2.  See  Appendix  C  for  explanation  of  abbreviations 


Searcher  Propagation  Loss  Curve  for  Case  A 
Examples . 
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Figure  3.3.  Searcher  Track  for  Case  A  Examples. 
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Searcher  Position  at  time  of  detection  for  A-I 


Figure  3.11.  Cumulative  Probability  of  Counter-Detection 
A-III. 
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C.  BARRIER  SEARCH:  SOME  FACTORS  AFFECTING  PENETRATION 

SUCCESS 

1.  Scenario  Description 

A  patrolling  submarine  (searcher)  has  set  up  a 
stationary  barrier  to  intercept  transiting  submarines 
(target) .  Frequently  submarines  exhibit  aspect  dependent 
radiated  noise  patterns  which  may  become  a  significant 
factor  in  a  barrier  scenario  where  the  searcher  is  most 
frequently  presented  a  bow-on  aspect.  It  is  desired  to  gain 
an  understanding  of  the  magnitude  of  this  problem,  and  to 
investigate  to  what  extent  the  distribution  of  targets 
across  the  barrier  will  affect  intercept  performance. 

2.  PASS  Input  Data 


Table  III  shows  the  input  data  for  the  base  case 


(B-I) . 


3.  Environment  and  Search  Plan 


Figures  3.18  and  3.19  show  the  propagation  loss 
curves  for  the  searcher  and  target.  Since  there  are  no 
convergence  zones,  the  MSEL  model  is  not  used. 

Figure  3.20  shows  the  search  area  to  be  a  90  nm  by 
200  nm  rectangle  with  the  searcher  conducting  a  stationary 
(back-and- forth)  barrier  perpendicular  to  the  expected 
direction  of  target  motion.  The  placement  of  the  barrier 
was  based  on  an  expected  median  detection  range  of  25  nm. 


4.  Platform  Parameters 

For  the  base  case  (B-I) : 

1.  The  searcher  speed  is  constant  at  10  knots. 

2.  The  searcher  FOH  is  87  db  and  is  not  target  aspect 
dependent . 

3.  The  initial  target  distribution  across  the  barrier  is 
uniform. 

4.  The  target  starts  at  the  top  of  the  barrier  (y=200) 
(i.e.  the  searcher  is  not  time  late  on  the  barrier) . 

5.  The  maximum  target  course  variation  is  plus-or-minus 
30  degrees. 

6.  The  target  speed  is  uniform  between  8  and  20  knots. 

7.  The  target  changes  course  and  speed  (independently) 
on  the  average  of  twice  per  hour. 

8.  The  target  FOM  is  84  db. 

5.  Results  and  Conclusions 

Numerical  results  for  the  base  case  (B-I)  and  varia¬ 
tions  of  the  base  case  (Case  B-II  and  B-III)  are  shown  in 
Table  IV.  The  variations  in  the  base  case  were: 

1.  Case  B-II:  The  distribution  of  the  initial  target 
position  across  the  top  of  the  barrier  was  altered  so 
that  50%  of  the  targets  were  uniformly  distributed 
between  0  and  30  miles  from  the  left  boundary,  and 
50%  of  the  targets  were  uniformly  distributed  between 
0  and  30  miles  from  the  right  boundary. 

2.  Case  B-III:  The  target  distribution  is  as  in  B-II. 
Additionally,  the  target  was  assumed  to  have  an 
aspect  dependent  source  level  which  induced  the 
target  aspect  dependent  searcher  FOM  as  shown  in 
Table  V.  The  significance  of  this  radiated  noise 
pattern  is  the  bow-null  presented  the  searcher  during 
barrier  approach  and  penetration. 
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The  simulation  results  indicate  that: 

1.  If  the  searcher  is  conducting  a  constant  speed  back- 
and-forth  barrier,  an  increase  in  the  probability  of 
penetration  can  be  realized  by  transiting  the  edge 
of  the  barrier. 

2.  When  calculating  search  effectiveness  in  a  barrier 
scenario,  the  aspect  dependency  of  target  source 
levels  may  be  a  significant  factor. 

A  complete  set  of  graphical  output  from  PASS  for 
Case  B-I  is  shown  in  Figures  3.21  through  3.25;  for  Case 
B— II  in  Figures  3.26  through  3.30;  for  Case  B-III  in 
Figures  3.31  through  3.39. 
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TABLE  IV 


NUMERICAL  RESULTS 

FOR  CASE 

B  EXAMPLES 

Case  B— I 

Case  B— II 

Case  B-III 

PD 

0.7055 

0.6040 

0.4410 

PDDP 

I. 0000 

1.0000 

1.0000 

MOES 

0.8521 

0.8014 

0.7152 

PCD 

0.0000 

0.0005 

0.0075 

PCDDP 

0.0000 

0.0000 

1.0000 

MOST 

0.2926 

0.3949 

0.5545 

ER 

2.91 

2.03 

1.29 

(T)S 

12.23 

12.39 

12.80 

(T)t 

none 

none 

11.73 

CR)3 

23.46 

24.25 

21.40 

<R-50)s 

22.99 

24.04 

21.05 

(R)  t 

none 

none 

5.83 

CR-50) t 

none 

none 

4.69 

Notes:  1. 

See  Appendix  C  for  explanation  of 

abbreviations . 

TABLE  V 

SEARCHER  FOM  VS.  TARGET 

ASPECT 

Target  Aspect 

Searcher  FOM 

000 

82 

045 

87 

090 

85 

135 

87 

180 

82 

225 

87 

270 

85 

315 

87 
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Searcher  Position  at  Time  of  Detection  for  B-I 


Figure  3.28.  Probability  of  Detection  Vs.  Range  for  B-II. 
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Figure  3.35. 


Searcher  Position  at  Time  of  Detection  for 
B— III. 
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Figure  3.38.  Target  Position  at  Time  of  Counter-Detection 
for  B-III . 
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D.  DEFINITE  RANGE  LAW  APPROXIMATIONS 

1.  Area  Search:  Best  Patrol  Speed  to  Evade  a  Systematic 

Searcher 

a.  Scenario  Description 

For  Case  C-I  the  searcher  is  given  a  definite 
range  lav  sensor  (detection  range  of  20  nn) ,  and  the 
stationary  target  (uniformly  distributed  over  the  search 
area.  A)  has  no  capability  to  counter-detect.  The  searcher 
systematically  sweeps  the  search  area  at  a  constant  speed 
(v)  with  track  spacing  equal  to  the  detection  range.  See 
Figure  3.24.  After  88  hours  of  search  (at  5  knots)  the 
searcher  has  come  within  one  detection  range  of  every  point 
in  the  search  area,  or  has  completely  swept  the  area. 

In  Case* I I  target  motion  is  introduced.  The 
target  travels  at  a  constant  speed  of  5  knots,  and  changes 
course  randomly  on  the  average  of  twice  an  hour. 

b.  PASS  Input  Data 

The  input  data  for  Case  C-I  is  shown  in  Table  VI 
Note  that  in  these  examples  figure-of-merit  is  chosen  to 
give  detection  ranges  based  on  propagation  loss  curves  shown 
in  Figures  3.40  and  3.41. 

c.  Results  and  Conclusions 

Numerical  results  for  Case  C-I  and  Case  C-II  are 


shown  in  Table  VII.  Graphical  results  for  Case  C-I  are 
shown  in  Figures  3.43  through  3.46,  and  for  Case  C-II  in 
Figures  3.47  through  3.49. 


A  comparison  of  the  results  shows: 

1.  The  times  to  detection  are  roughly  uniform  for  Case 
C-I  as  can  be  seen  in  the  very  nearly  constant  slope 
of  Figure  3.43.  Note  the  departure  from  exhaustive 
search  results  (manifested  by  the  jump  at  time  *  0  in 
Figure  3.43],  which  is  due  to  the  significant  prob¬ 
ability  of  detection  at  the  instant  the  search 
starts.  Also,  "bumps"  in  the  cumulative  probability 
of  detection  curve  occur  at  turn  points,  and  repre¬ 
sent  the  combined  effects  of  the  searcher  having  a 
significant  portion  of  his  detection  circle  outside 
the  search,  area,  and  sweeping  area  previously  swept, 
thus  resulting  in  a  low  detection  rate. 

2.  Introducing  target  motion  in  Case  C-II  skews  the 
distribution  of  time  to  detection  to  the  right  (see 
Figure  3.47],  resulting  in  am  increased  mean  time 
to  detection. 

Conclusion  2  above  raises  the  question:  What  is 
the  "best"  speed  to  evade?  The  results  of  Case  C-II  indicate 
that  the  best  speed  is  not  zero,  if  the  criteria  is  mean 
time  to  detection.  Before  proceeding,  the  following  terms 
are  defined: 

1.  Speed  Ratio:  The  ratio  of  target  speed  (u)  to 
searcher  speed  (v] .  Note  that  if  the  speed  ratio  is 
zero,  then  the  target  is  stationary. 

2.  Mean  Time  Ratio:  Hie  ratio  of  the  mean  time  to  detec¬ 
tion  (for  a  given  speed  ratio)  to  the  mean  time  of 
detection  when  the  speed  ratio  is  zero. 

3.  Median  Time  Ratio:  Hie  ratio  of  the  median  time  to 
detection  (for  a  given  speed  ratio)  to  the  median 
time  to  detection  when  the  speed  ratio  is  zero. 

To  investigate  the  problem  of  best  evasion 
speed,  PASS  was  run  for  various  search  speeds  (  v  =  5,  10, 

15,  20  knots)  and  various  speed  ratios  at  each  search  speed, 
under  case  C  conditions  (definite  range  law  sensor) . 


Figures  3.50  and  3.51  show  the  results  of  these 
simulations.  Note  that  whenever  the  time  ratio  is  greater 
them  1.0,  the  time  to  detection  is  greater  them  the  time  to 
detection  for  a  stationary  target.  Figure  3.50  shows  the 
meem  time  to  detection  peaks  at  speed  ratios  between  1.0  and 
1.4,  depending  on  the  searcher  speed.  Figure  3.51  shows  the 
mediem  time  to  detection  is  relatively  insensitive  to  speed 
ratios  below  1.0,  and  drops  off  rather  sharply  for  speed 
ratios  above  l.Q.  The  conclusion  drawn  is  that  the  best 
speed  to  evade  the  searcher  is  about  equal  to  the  searcher 
speed.  If  the  searcher  speed  is  not  known,  the  best  tactic 
would  be  to  limit  maximum  patrol  speed  to  the  most  likely 
minimum  searcher  speed.  It  is  important  to  note  that  the 
best  patrol  speed  is  not  necessarily  the  minimum  achievable 
speed. 
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TABLE  VII 

NUMERICAL  RESULTS  FOR  CASES  C-I  AND  C-II 


Case  C-I 

Case  C-II 

PD 

1.0000 

0.9992 

PDDP 

1.0000 

1.0000 

HOES 

1.0000 

0.9994 

PCD 

0.0000 

0.0000 

PCDDP 

0.0000 

0.0000 

MOET 

0.0000 

0.0002 

.  ER 

infinite 

4496.00 

(T)S 

40.12 

43.84 

Var(T)s 

707.17 

1195.38 

TCnin) 

0.00 

0.00 

TCnax) 

87.24 

302.07 

Notes:  1.  See  Appendix  C  for  explanation  of 
abbreviations . 


a. 


Figure  3.40.  Searcher  Propagation  Loss  for  Case  C  Examples 


Figure 
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.41.  Target  Propagation  Loss  for  Case  C  Examples. 
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E.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 

In  the  proceeding  examples,  it  was  intended  to  demon¬ 
strate  the  capabilities  and  flexibility  of  the  PASS  program. 
Of  significance  is  the  fact  that  an  experienced  operator, 
using  the  interactive  data  entry  procedures,  could  run  all 
of  the  example  problems  in  about  four  hours  (including  time 
for  hard-copy  printing  and  graphics  production) .  These  runs 
represent  274,000  individual  simulation  replications. 

The  utility  of  the  graphical  output  should  be  evident. 
Frequently  the  pictorial  displays  provide  insight  that  is 
obscured  by  a  large  quantity  of  numbers.  Additionally,  the 
very  shape  of  certain  graphical  results  inspire  questions 
the  answer  to  which  are  critical  to  a  thorough  analysis. 
Improvements  or  additions  to  PASS  might  include: 

1.  Multiple  sensor  capability  by  both  platforms. 

2.  Multiple  search  platforms  (e.g.  searcher,  tax  get, 
and  consort) . 

3.  Post  detection  analysis  (e.g.  break  contact, 
approach  and  attack,  target  motion  analysis) . 

4.  Gaming  capability  (real  time  control  over  platform 
motion) . 

5.  Weapon  employment  capability. 


APPENDIX  A 


PASS  USERS  GUIDE 


A.  INTRODUCTION 

This  User's  Guide  is  a  set  of  instructions  designed  to 
allow  the  execution  of  program  PASS  on  the  NPGS  computer  by 
users  unfamiliar  with  the  program  code.  The  user  should  be 
familiar  with  the  concepts  of  the  simulation  model  as 
presented  in  Chapter  I,  and  have  available  for  reference  the 

variable  list  in  Appendix  B. 

» 

1.  Data  File  Management 

When  PASS  is  run,  an  input  data  file  may,  or  may  not 
exist.  If  no  data  file  exists,  the  easiest  way  to  create  one 
is  to  execute  program  PASS  and  input  data  during  the  inter¬ 
active  session.  If  a  data  file  exists  it  must  have  the 
filename  PASS  and  filetype  DATA.  If  a  data  file  is  inter¬ 
actively  created,  it  will  be  in  file  PASS  DATA  at  the 
completion  of  the  run  of  PASS. 

At  the  start  of  each  run,  if  PASS  DATA  exists,  it  is 

copied  into  file  PASS  HOLD.  If  changes  are  made  to  the  input 

data,  file  PASS  DATA  will  reflect  these  changes,  and  PASS 

HOLD  will  contain  the  unchanged  data.  If  no  changes  were 

made,  PASS  HOLD  and  PASS  DATA  will  be  identical.  To  save  the 

data  in  PASS  HOLD,  it  must  be  renamed  before  the  next  run  of 

PASS.  PASS  DATA  will  remain  unchanged  as  long  as  no  changes 

are  made  to  the  data  during  the  interactive  session. 
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Upon  completion  of  each  PASS  run  a  file  PASS  OUTPUT 
is  created.  Depending  on  the  option  selected  by  the  user, 
this  file  will  contain  some,  or  all  of  the  following: 

1.  A  complete  record  of  the  input  data,  or  a  user 
specified  run  identification  number. 

2.  A  summary  of  the  results  of  the  simulation. 

3.  Sectioning  results  of  times  and  ranges  of  detection 
and  counter-detection  (requires  a  minimum  of  500 
data  points) . 

4.  Histogram  and  statistics  of  times  and  ranges  of 
detection  and  counter-detection  (requires  a  minumum 
of  10  data  points) . 

If  this  file  is  to  be  saved,  it  must  be  renamed  or  printed 
before  the  next  run  of  PASS. 

Upon  completion  of  each  run  of  PASS  a  file  PASPLT 
DATA  is  created.  This  file  provides  the  data,  in  the  proper 
format,  to  be  used  by  the  FORTRAN  program  PASPLT.  Program 
PASPLT  uses  the  DISSPLA  graphics  system  resident  on  the  NPGS 
computer  to  produce  a  graphical  display  of  the  simulation 
input  parameters  and  results. 

2 .  Required  Input  Data 

PASS  requires  the  following  data,  which  the  user 
should  have  ready  for  interactive  data  input  if  a  complete 
PASS  DATA  file  does  not  exist: 

1.  Size  of  the  search  area  rectangle. 

2.  Searcher  track  anchor  points  and  return  point. 

3.  Searcher  propagation  loss  data  (direct  path  and  CZ's). 

4.  Target  propagation  loss  data  (direct  path  and  CZ's). 


5.  Target  initial  x-distnbution  if  not  uniform  (barrier 
scenario  only) . 

6.  Searcher  figure-of-merit  at  sprint  and  drift  speeds 
(as  a  function  of  target  aspect,  if  target  source 
level  is  aspect  dependent) . 

7.  Searcher  sprint  and  drift  speeds. 

8.  Searcher  sprint  and  drift  times. 

9.  Target  figure-of-merit  against  sprinting  and  drifting 
searcher . 

10.  Target  speed  range. 

11.  Rate  of  target  zigs  (course  and  speed  changes) . 

12 .  Acousti  '  fluctuation  rate  parameters . 

13.  Acoustic  fluctuation  scale  parameters. 

14.  Run  parameters  (random  number  seed,  number  of  repli¬ 
cations  ,  maximum  search  time  for  each  replication). . 


During  the  interactive  input  of  data,  the  program 
does  some  error  checking.  However,  it  is  far  from  impossible 
to  enter  bad  data',  or  incorrectly  formatted  data.  If  the 
program  fails  to  run  as  expected,  the  most  probable  cause  is 
bad  data. 

3.  The  Executive  Program 

The  executive  program  PASS  EXEC  is  necessary  to  run 
the  FORTRAN  program  PASS.  The  executive  program  conducts 
data  file  management  and  starts  PASS.  Prior  to  invoking  the 
executive  program,  PASS  must  be  compiled  in  either  FORTRAN-G 
or  FORTRAN-H  with  the  text  file  on  the  A-DISK.  The 
FORTRAN-H  compiler  is  recommended  as  it  results  in  a  faster 
program  run  time. 


B.  PROGRAM  PASS  EXECUTION 


Each  of  the  following  sections  describes  a  terminal 
option.  .  Each  option  presents  a  new  menu  on  the  terminal 
screen.  Depending  on  your  response  to  each  option  menu,  you 
will  proceed  to  the  next  option,  or  you  will  be  prompted  to 
enter  data  or  make  some  decision  regarding  the  nature  of  the 
simulation. 

1.  Starting  the  Program:  Invoking  the  Executive 

Program 

To  start  program  PASS,  ENTER:  PASS 
The  screen  will  contain  a  quick  review  of  data-file  manage¬ 
ment  and  the  first  OPTION. 

2.  OPTION  NO.  1:  Terminate  or  Proceed 

This  option  allows  you  to  "gracefully"  terminate 
the  program. 

To  terminate  the  program,  ENTER;  1 

To  proceed  with  execution,  ENTER:  2 
If  you  terminate  now,  and  you  had  an  input  data  file  PASS 
DATA,  it  will  have  been  copied  into  file  PASS  HOLD,  and  PASS 
DATA  will  be  empty. 

3.  MASTER  OPTION:  Accept  Program  Defaults 

This  option  allows  you  to  bypass  the  interactive 
sections  and  proceed  directly  to  the  simulation.  Bypassing 
requires  a  complete  and  properly  formatted  input  data  file, 
and  the  acceptance  of  the  program  defaults. 


SjMXAHSf  iS&KXfiCa&L  StettSfe 


--4  ’>  V  ■>  -.4  V  "/  V 


To  list  the  program  defaults,  ENTER:  1 
To  accept  the  program  defaults,  ENTER:  2 
To  commence  interactive  options,  ENTER:  3 

4.  OPTION  NO.  2A:  Signal  Integration  Model 
This  option  allows  you  to  select  a  threshold 

crossing  model  without  signal  integration,  or  the  MSEL 
integration  model.  It  is  recommended  that  the  MSEL  model 
be  used  if  strong  convergence  zones  are  present. 

To  select  no  signal  integration,  ENTER:  1 
To  select  the  MSEL  integration  model,  ENTER:  2 

5.  OPTION  NO.  2:  Acoustic  Fluctuation  Model 
This  option  allows  you  to  select  either  a 

Lambda-Sigma  Jump  or  a  Gauss-Markov  error  process  to  model 
acoustic  fluctuations. 

To  select  the  Lambda-Sigma  Jump  model,  ENTER:  1 
To  select  the  Gauss-Markov  model,  ENTER:  2 

6.  OPTION  NO.  3A:  Existence  of  File  PASS  DATA 

The  program  must  be  "told"  whether  or  not  to  read 
data  from  an  existing  data  file. 

If  PASS  DATA  is  not  on  your  A-DISK,  ENTER:  1 
If  PASS  DATA  is  on  your  A-DISK,  ENTER:  2 

7.  OPTION  NO.  3:  Search  Area  Size 
Variables:  XMAX,  YMAX 

This  option  allows  you  to  specify  the  size  of  the 
search  area  in  which  the  target  is  confined.  The  area  dimen 
sions  are  in  nautical  miles. 


VST 


To  review  existing  data,  ENTER:  1 


To  accept  existing  data,  ENTER:  2 


To  change  existing  data,  ENTER:  3 


If  (1)  is  entered,  the  existing  data  (in  file  PASS 


DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 


this  data,  ENTER:  0. 


If  (2)  is  entered,  the  existing  data  will  be  read 


into  PASS,  and  the  next  option  menu  will  be  presented  on 


the  screen. 


If  (3)  is  entered,  prompting  messages  will  be  on  the 


terminal  screen  to  allow  you  to  enter  new  data.  When  data 


entry  is  completed,  you  must  return  to  the  option  menu  (to 


review,  accept,  or  further  modify  data).  To  return. 


ENTER:  0 . 


8.  OPTION  NO.  4:  Searcher  Track 


Variables:  NP,  KP,  XP(I) ,  YP(I) 


This  option  allows  you  to  enter  the  searcher  track 


anchor  points  (2  to  50  points) ,  and  the  return  point.  The 


dimensions  of  the  anchor  points  are  nautical  miles  from 


(0,  0) 


To  review  existing  data,  ENTER:  1 


To  accept  existing  data,  ENTER:  2 
To  change  existing  data,  ENTER:  3 


If  (1)  is  entered,  the  existing  data  (in  file  PASS 


DATA  on  your  A-DISK)  will  be  presented  on  the  screen 
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for  review.  To  return  to  the  option  menu  to  accept  or 
modify  this  data,  ENTER:  0. 


If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  messages  will  be  on  the 
terminal  screen  to  allow  you  to  enter  new  data .  When  data 
entry  is  completed,  you  must  return  to  the  option  menu  (to 
review,,  accept,  or  futher  modify  data)  .  To  return, 

ENTER:  0. 

9.  OPTION  NO.  5:  Searcher  Direct  Path  Propagation  Loss 

Variables:  MO,  RO(I) ,  0L(I) 

This  option  allows  you  to  input  the  propagation  loss 
of  the  searcher  as  a  function  of  range.  The  number  of  points 
should  be  between  2  and  20,  ensuring  that  there  are  suffi¬ 
cient  points  to  adequately  describe  the  propagation  loss 
curve.  Range  is  in  nautical  miles,  and  propagation  loss  in 
decibels  (db) . 

To  review  existing  data,  ENTER:  1 

To  accept  existing  data,  ENTER:  2 

To  change  existing  data,  ENTER:  3 

If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 


If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  messages  will  be  on  the 
terminal  screen  to  allow  you  to  enter  new  data.  When  data 
entry  is  completed,  you  must  return  to  the  option  menu  (to 
review,  accept,  or  further  modify  data) .  To  return, 

ENTER:  0. 

10.  OPTION  NO.  6:  Target  Direct  Path  Propagation  Loss 

Variables:  MT,  RT(I)  ,  OL(I) 

This  option  is  similar  to  OPTION  5,  but  here  you 
enter  the  propagation  loss  data  for  the  target. 

To  review  existing  data,  ENTER:  1 

To  accept  existing  data,  ENTER:  2 

To  change  existing  data,  ENTER:  3 

If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  messages  will  be  on  the 
terminal  screen  to  allow  you  to  enter  new  data.  When  data 
entry  is  completed,  you  must  return  to  the  option  menu  (to 
review,  accept,  or  further  modify  data) .  To  return,  ENTER:  0 
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11.  OPTION  NO.  7;  Searcher  Convergence  Zones 
Variables:  NCZOf  RCZO(I,J) ,  CZLO(I) 

This  option  allows  you  to  input  the  number  of 
searcher  convergence  zones  (0  to  5) ,  the  width  of  each  CZ, 
and  the  propagation  loss  in  each  CZ.  The  CZs  are  modeled  as 
inverted  square-wells  superimposed  on  the  OP  propagation 
loss  curve.  You  will  have  to  make  a  subjective  interpreta¬ 
tion  of  the  actual  CZ  data  to  determine  the  range  to  the 
inner  and  outer  rings  of  the  square-well,  and  the  effective 
propagation  loss  therein.  Ranges  are  in  nautical  miles, 
and  propagation  loss  in  db. 

To  review  existing  data,  ENTER:  1 
To  accept  existing  data,  ENTER:  2 
To  change  existing  data,  ENTER:  3 

If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  messages  will  be  on 
the  terminal  screen  to  allow  you  to  enter  new  data.  When 
data  entry  is  completed,  you  must  return  to  the  option  menu 
(to  review,  accept,  or  further  modify  data) .  To  return. 


ENTER:  0 . 


Variables:  NCZT,  RCZT(I,J) ,  CZLT(I) 

This  option  is  similar  to  OPTION  7,  but  here  you 
enter  the  data  for  target  convergence  zones. 

To  review  existing  data,  ENTER:  1 
To  accept  existing  data,  ENTER:  2 
To  change  existing  data,  ENTER:  3 
If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  messages  will  be  on 
the  terminal  screen  to  allow  you  to  enter  new  data.  When 
data  entry  is  completed,  you  must  return  to  the  option 
menu  (to  review,  accept,  or  further  modify  data) .  To 
return,  ENTER:  0. 

13.  OPTION  NO.  9:  Area  or  Barrier  Search 

This  option  allows  you  to  select  either  an  area  or 
barrier  search  scenario.  If  area  search  is  selected,  the 
next  option  will  be  10A.  If  barrier  search  is  selected,  the 
next  option  will  be  9A. 

To  select  area  search  scenario,  ENTER:  1 
To  select  barrier  search  scenario,  ENTER:  2 
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14.  OPTION  NO.  9A:  Initial  Target  Lateral  Distribution 

Variables:  NBINS,  XBIN(I) ,  PBIN(I) 

This  option  allows  you  to  specify  the  lateral 
distribution  of  target  starting  position  (barrier  scenario 
only) .  Distances  are  in  nautical  miles  from  X  *  0 . 

To  distribute  the  initial  positions  uniformly  on 
(0 , XMAX) ,  ENTER:  1. 

To  specify  the  lateral  distribution,  ENTER:  2 

15.  OPTION  NO.  10A:  Searcher  FOM  Dependent  on  Target 

Aspect 

Variables:  NLS ,  BRG(I) ,  FOMBD(I) ,  FOMBS(I) 

This  option  allows  you  to  specify  the  searcher 
figure-of-merit  (at  drift  and  sprint  speeds)  as  a  function 
of  the  relative  bearing  from  target  to  searcher.  This  is 
equivalent  to  a  target  source  level  dependent  upon  aspect 
angle.  Angles  are  in  degrees,  and  FOM  in  db. 

To  specify  searcher  FOM  not  target  aspect  depen¬ 
dent,  ENTER:  1 

To  specify  searcher  FOM  target  aspect  dependent, 

ENTER:  2 

If  (1)  is  selected,  OPTION  10  will  be  presented  next. 


If  (2)  is  selected,  after  data  entry,  OPTION  11  will  be 
presented  next. 


16.  OPTION  NO.  10:  Searcher  Figure-of-Merit 
Variables:  FOMOS ,  FOMOD 

This  option  allows  you  to  specify  the  searcher 
figure-of-merit  against  the  target  when  the  searcher  is  at 
sprint  or  drift  speeds.  FOM  IN  db. 

To  review  existing  data,  ENTER:  1 
To  accept  existing  data,  ENTER:  2 
To  change  existing  data,  ENTER:  3 
If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  messages  will  be  on 
the  terminal  screen  to  allow  you  to  enter  new  data.  When 
data  entry  is  completed,  you  must  return  to  the  option 
menu  (to  review,  accept,  or  further  modify  data) .  To 
return,  ENTER:  0. 

17.  OPTION  NO.  11:  Searcher  Sprint  and  Drift  Speeds 
Variables:  SOS,  SOD 

This  option  allows  you  to  select  the  search  speed 


(sprint  and  drift) .  For  constant  speed  search,  set  SOS  - 
SOD.  Speeds  in  knots. 


To  review  existing  data,  ENTER:  1 
To  accept  existing  data,  ENTER:  2 
To  change  existing  data,  ENTER:  3 
If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  menu  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  mess^ge^will  be  on 
the  terminal  screen  to  allow  you  to  enter  new  data.  When 
data  entry  is  completed,  you  must  return  to  the  option 
menu  (to  review,  accept,  or  further  modify  data) .  To 
return,  ENTER:  0. 

18.  OPTION  NO.  12:  Searcher  Sprint  and  Drift  Times 
Variables:  TS,  TD 

This  option  allows  you  to  specify  the  time  the 
searcher  spends  at  sprint  and  drift  speeds.  Time  in  hours. 
To  review  existing  data,  ENTER:  1 
To  accept  existing  data,  ENTER:  2 
To  change  existing  data,  ENTER:  3 
If  (1)  is  entered,  the  existing  data  (in  file  PASS 


DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 


If  (2)  is  entered ,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on  the 
screen. 

If  (3)  is  entered,  prompting  messages  will  be  on  the 
terminal  screen  to  allow  you  to  enter  new  data.  When  data 
entry  is  completed,  you  must  return  to  the  option  menu  (to 
review,  accept,  or  further  modify  data) .  To  return, 

ENTER:  0 . 

19.  OPTION  NO.  13:  Target  Figure-of-Merit 

Variables :  FOMTS ,  FOMTD 

This  option  allows  you  to  specify  the  target 
figure-of-roerit  against  the  searcher  when  the  searcher  is  at 
sprint  or  drift  speed.  FOM  in  db. 

To  review  existing  data,  ENTER:  1 

To  accept  existing  data,  ENTER:  2 

To  change  existing  data,  ENTER:  3 

If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on  the 
screen. 

If  (3)  is  entered,  prompting  messages  will  be  on 
the  terminal  screen  to  allow  you  to  enter  new  data.  When 
data  entry  is  completed,  you  must  return  to  the  option 
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menu  (to  review,  accept,  or  further  modify  data) .  To 
return,  ENTER:  0. 

20 .  OPTION  NO.  14:  Target  Speed  Range 

Variables:  STMIN,  STMAX 

This  option  allows  you  to  specify  the  range  of 
target  spped.  Target  speed  is  distributed  uniform  on 
(STMIN, STMAX) .  Speed  in  knots. 

To  review  existing  data,  ENTER:  1 

To  accept  existing  data,  ENTER:  2 

To  change  existing  data,  ENTER:  3 

If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  messages  will  be  on 
the  terminal  screen  to  allow  you  to  enter  new  data.  When 
data  entry  is  completed,  you  must  return  to  the  option  menu 
(to  review,  accept,  or  further  modify  data) .  To  return, 
ENTER:  0. 

21.  OPTION  NO.  15:  Rate  of  Target  Zigs 

Variables:  RTCC,  RTSC 

This  option  allows  you  to  specify  the  rate  of  target 


course  and  speed  changes.  Rates  are  "per-hour" . 


To  review  existing  data,  ENTER:  1 
To  accept  existing  data,  ENTER:  2 
To  change  existing  data,  ENTER:  3 

If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
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review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on  the 
screen. 

If  (3)  is  entered,  prompting  messages  will  be  on  the 
terminal  screen  to  allow  you  to  enter  new  data.  When  data 
entry  is  completed,  you  must  return  to  the  option  menu  (to 
review,  accept,  or  further  modify  data) .  To  return, 

ENTER:  0. 

22.  OPTION  NO.  16:  Acoustic  Fluctuation  Rate  Parameters 
Variables:  ALAM(I) 

This  option  allows  you  to  specify  the  rate  param¬ 
eters  for  the  acoustic  fluctuation  process.  Rates  are 
"per-hour" . 

To  review  existing  data,  ENTER:  1 
To  accept  existing  data,  ENTER:  2 
To  change  existing  data,  ENTER:  3 
If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen 


for  review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 


If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  messages  will  be  on 
the  terminal  screen  to  allow  you  to  enter  new  data.  When 
data  entry  is  completed,  you  must  return  to  the  option  menu 
(to  review,  accept,  or  further  modify  data) .  To  return, 
ENTER:  0. 

23.  OPTION  NO.  17:  Acoustic  Fluctuation  Scale 

Parameters 

Variables:  SIGMA (I) 

This  option  allows  you  to  specify  the  scale  param¬ 
eters  for  the  acoustic  fluctuation  process.  Scale  in  db. 

To  review  existing  data,  ENTER:  1 

To  accept  existing  data,  ENTER:  2 

To  change  existing  data,  ENTER:  3 

If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 


the  screen. 


If  (3)  is  entered,  prompting  messages  will  be  on 


the  terminal  screen  to  allow  you  to  enter  new  data.  When 
data  entry  is  completed,  you  must  return  to  the  option  menu 
(to  review,  accept,  or  futher  modify  data) .  To  return, 

ENTER:  0. 

24 .  OPTION  NO.  18:  Run  parameters 

Variables:  SEED,  NREP,  TMAX 

This  option  allows  you  to  select  parameters  that 
govern  the  simulation  execution  and  random  number  generation. 
Time  in  hours. 

To  review  existing  data,  ENTER:  1 

To  accept  existing  data,  ENTER:  2 

To  change  existing  data,  ENTER:  3 

If  (1)  is  entered,  the  existing  data  (in  file  PASS 
DATA  on  your  A-DISK)  will  be  presented  on  the  screen  for 
review.  To  return  to  the  option  menu  to  accept  or  modify 
this  data,  ENTER:  0. 

If  (2)  is  entered,  the  existing  data  will  be  read 
into  PASS,  and  the  next  option  menu  will  be  presented  on 
the  screen. 

If  (3)  is  entered,  prompting  messages  will  be  on 
the  terminal  screen  to  allow  you  to  enter  new  data.  When 
data  entry  is  completed,  you  must  return  to  the  option 
menu  (to  review,  accept,  or  further  modify  data) .  To 
return,  ENTER:  0. 
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25 .  OPTION  NO.  19:  Input  Data  to  Output  File 
Hiis  option  allows  you  to  either  send  all  of  the 

input  data  to  the  output  file,  or  to  identify  the  output 
file  with  a  run  identification  number. 

To  send  all  input  data  to  the  output  file, 

ENTER:  1 

To  supress  input  data  echo  and  select  a  run  ID, 

ENTER:  2 

26 .  OPTION  NO.  20:  Statistical  Analysis 
This  option  allows  you  to  select  the  amount  and 

type  of  statistical  analysis  to  be  done  on  the  simulation 
results . 

For  complete  statistical  analysis  (sectioning, 
histograms,  and  parameter  point  estimates,  ENTER:  1 
To  omit  sectioning,  ENTER:  2 
For  no  statistical  analysis,  ENTER:  3 

C.  NORMAL  RUNNING  INDICATION 

Immediately  following  OPTION  20,  the  screen  will 
display  a  running  total  of  the  number  of  replications 
completed.  This  display  is  updated  every  200  replication, 
which  allows  the  user  to  follow  the  progress  of  program 
execution.  Upon  completion  of  the  last  replication,  the 
screen  will  be  cleared,  and  the  simulation  quick-summary 
will  be  presented. 
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D.  GRAPHICAL  PRESENTATION  OF  RESULTS:  PROGRAM  PASPLT 
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The  FORTRAN  program  PASPLT  is  intended  to  be  used  in 
conjunction  with  the  DISSPLA  graphics  system  to  produce 
graphical  representations  of  the  results  of  a  run  of  program 
PASS.  To  run  program  PASPLT  you  will  need: 

1.  Data  file  PASPLT  DATA  on  your  A-DISK.  PASPLT  DATA  is 
automatically  created  by  running  program  PASS . 

2.  Dual  terminal  graphics  capability  (TEKTRONIX  618  and 
thermal  printer) . 

3.  Program  PASPLT  FORTRAN-G  or  FORTRAN-H  compiled  on 
your  A-DISK. 

To  produce  graphical  output,  invoke  the  graphics  package 
by  entering:  DISSPLA.  Two  responses  to  the  DISSPLA  execu¬ 
tive  program  prompts  are  required: 

1.  You  must  identify  the  program  PASPLT  as  the  FORTRAN 
program  to  be  used  with  DISSPLA. 

2.  You  must  define  the  data  file  as  follows:  FILEDEF  10 
DISK  PASPLT  DATA. 

The  above  is  done  in  response  to  prompting  messages  produced 
on  the  screen  after  DISSPLA  is  entered. 

Each  plot  is  presented  on  the  TEKTRONIX  618  screen,  and 
hard  copies  may  be  made  by  pressing  the  HARD  COPY  button  at 
the  bottom  of  the  screen.  To  produce  the  next  plot,  hit 
ENTER  on  the  keyboard.  If  any  messages  (other  them  error 
messages)  appear  on  the  IBM  screen,  clear  the  screen  until 
the  next  plot  appears  on  the  TEKTRONIX  screen. 
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APPENDIX  B 


DESCRIPTION  OF  VARIABLE  NAMES  USED  IN  PASS 


The  following  list  of  major  variable  names  used  in  PASS 
is  arranged  alphabetically .  Those  variable  names  with  an 
asterik  (*)  indicate  that  an  understanding  of  their  function 
is  required  to  run  PASS  in  conjunction  with  Appendix:  A  (PASS 
User's  Guide).  Following  each  variable  name  is  the  data  type 
and  physical  dimension,  if  any,  and  then  a  verbal  description 
of  the  variable  use.  For  input  data,  restrictions  (e.g. 
variable  size,  entry  procedures)  are  listed  under  NRESTR:n. 
Some  variables  used  for  temporary  storage,  loop  counters, 
or  pointers  are  not  listed  as  their  function  should  be 
apparent  to  the  reader  of  the  code. 

PASS  VARIABLE  NAMES: 

1.  AFLtI) /real/decibels/. ... .Value  of  the  Ith  acoustic 
fluctuation  process.  1*1,2, 3. 

2.  AFLGM/logical/ . If  the  Gauss-Markov  error  process 

is  to  be  used ,  this  variable  is  TRUE .  Otherwise ,  it 
is  false. 

3.  AFLLS/logical/ . If  the  Lambda-Sigma  Jump  error 

process  is  used,  this  variable  is  TRUE.  Otherwise, 
it  is  FALSE. 

4.  ALAM(I) /*/real/inverse  hours/ . Rate  parameter  for 

the  Ith  fluctuation  process,  1*1, 2, 3,  (input). 

RESTR:  positive. 

5.  ALFA/ real/ . Multiplier  for  the  GMA  process. 


ANGl/real/degrees (external) ,  radians (internal) / . A 

limit  of  target  courses  allowed  in  the  barrier  scenario 
In  the  simulation,  courses  are  measured  CCW  from  the 
X-axis.  In  the  barrier  scenario,  the  target  base 
course  is  in  the  negative  Y-direction,  or  psuedo- 
course  of  270-degrees.  ANG1  *  270  -  ANG2. 

ANG2/*/real/degrees  (external)  r  radians  ( internal)  / . 

Angular  variation  from  base  course  (ANG1)  allowed  for 
the  target  in  the  barrier  scenario.  The  target  courses 
are  distributed  uniformly  over  psuedo-course  270±ANG2. 
See  ANG1  (input) .  RESTR:  non-negative,  not  greater 
them  90-degrees. 

ASPECT/logical/ . If  the  target  source  level  is 

aspect  dependent,  this  variable  is  TRUE.  Otherwise, 
it  is  FALSE. 

BETA/real/ . Multiplier  for  the  GMA  process. 

BREL/real/radians/ . Relative  bearing  from  the 

target  to  the  searcher. 

BRG(I)  /*/real/degrees  (external) ,  radians  (internal)  / 

. Relative  bearing  associated  with  target  aspect 

related  searcher  figure-of-merit,  I  *  1,...,NLS, 
(input).  RESTR:  BRG(l)  must  be  000.0  (directly  on 
the  bow)  ,  and  all  subsequent  BRG(I)  must  be  in 
ascending  order,  such  that  BRG(I)  <  BRG(I+1) ,  and 
BRG(I)  max  *  360.0. 

Cl,C2,C3,C4/integer/ . Counters  for  STACK1,2 , 3,4. 

COSX/real/ . X-component  of  a  unit  vector  in  the 

direction  of  target  motion. 

CZLO(I) /*/real/decibels/ . Propagation  loss  in  the 

Ith  searcher  convergence  zone  (input) .  If  the  number 
of  searcher  convergence  zones  (NCZO)  is  zero,  then 
this  data  is  not  input.  Otherwise,  1*1, . . . ,NCZO. 

RESTR :  pos i t i ve . 

CZLT (I) /*/real/decib61s/ . Propagation  loss  in  the 

Ith  target  convergence  zone  (input) .  If  the  number 
of  target  convergence  zones  (NCZT)  is  zero,  then  this 
data  in  not  input.  Otherwise,  1*1,... NCZT. 

RESTR:  positive. 

Dl/real/nautical  miles/ . Distance  from  current 

searcher  position  (XO,YO)  to  the  projected  searcher 
position  (XOT, YOT) . 


17.  D2/real/nautical  miles/ . Distance  from  current 

searcher  position  (XO,YO)  to  the  end  of  the  current 
search  leg. 

18.  DIST( I) /real/nautical  miles/ . The  distance  from 

searcher  track  anchor  point  I  to  anchor  point  (1+1) . 
The  length  of  the  Ith  leg  of  the  search  pattern. 

19.  DX(I)/real/ . X-component  of  the  unit  vector  in  the 

direction  of  searcher  motion  from  anchor  point  I  to 
anchor  point  (.1+1)  . 

20.  DY(I)/real/ . Y-component  of  the  unit  vector 

pointing  in  the  direction  of  searcher  motion  from 
anchor  point  I  to  anchor  point  (1+1) . 

21.  ECHOS/logical/ . If  the  user  chooses  to  have  the 

input  data  sent  to  the  output  file,  this  variable  is 
TRUE.  Otherwise,  it  is  FALSE. 

22.  ENUF/logical/ . If  5  hours  of  figure-of-merit  data 

for  the  target  has  been  recorded,  this  variable  is 
TRUE.  Otherwise,  it  is  FALSE. 

23.  ENUFF/logical/ . If  5  hours  of  figure-of-merit  data 

for  the  searcher  has  been  recorded,  this  variable  is 
TRUE.  Otherwise,  it  is  FALSE. 

24.  FO(I) /real/decibels/ . A  sample  of  the  searcher 

figure-of-merit  for  the  first  five  hours  of  a  repli¬ 
cation  used  to  provide  a  representative  plot  of  FOM 
as  a  reference.  This  variable  is  not  used  computa¬ 
tionally  in  the  program. 

25.  FOMBD(I) /*/real/decibels/ . Searcher  figure-of- 

merit  when  at  drift  speed  and  with  a  relative  bearing 
from  the  target  to  the  searcher  of  BRG(I) ,  where 

I  =  1 , . . . , NLS ,  (input).  RESTR:  positive. 

26 .  FOMBS (I) /*/real/decibels/ . Searcher  figure-of- 

merit  at  sprint  speed  at  the  corresponding  relative 
bearing  from  target  to  searcher  of  BRG(I) ,  I  =  1,..., 
NLS,  (input) .  RESTR:  positive. 

27.  FOMO/real/decibels/ . Current  searcher  figure-of- 

merit. 

28.  FOMOD/*/real/decibels/ . The  searcher  figure-of- 

merit  at  drift  speed  (input) .  RESTR:  positive. 


141 


29.  FQMOS/*/real/decibels/ . The  searcher  figure-of- 

merit  at  sprint  speed  (input) .  RESTR:  positive. 

30.  FOHT/real/decibels/ . Current  target  figure-of- 

merit. 

31.  FOMTD/*/real/decibels/ . The  target  figure-of-merit 

when  the  searcher  is  at  drift  speed  (input) . 

RESTR:  positive. 

32.  FOMTS/*/real/decibels/ . The  target  figure-of-merit 

when  searcher  is  at  sprint  speed  (input) . 

RESTR:  positive. 

33.  FILE/logical/ . This  variable  is  TRUE  if  the  user 

specifies  a  data  file,  PASS  DATA,  presently  exists  on 
his  A-disk.  Otherwise,  it  is  FALSE. 

34.  FT (I) /real/decibels/ . A  sample  of  the  target 

figure-of-merit  for  the  first  five  hours  of  a  repli¬ 
cation  used  to  provide  a  plot  of  representative  FOM 
as  a  reference.  ‘  This  variable  is  not  used  computa¬ 
tionally  elsewhere  in  the  program. 

35.  G04IT/Iogical/ . This  variable  is  TRUE  if  the  user 

specifies  that  the  interactive  data  input  is  to  be 
bypassed  and  input  data  is  to  be  read  directly  from 
an  existing  file  on  A-disk,  PASS  DATA.  Otherwise, 
it  is  FALSE. 

36.  ICZO/integer/ . Identifies  the  number  of  the 

convergence  zone  in  which  the  searcher  CZ  detection 
took  place. 

37.  ICZT/integer/ . Identifies  the  number  of  the 

convergence  zone  in  which  the  target  CZ  counter¬ 
detection  took  place. 

38.  INTEG/logical/ . This  variable  is  TRUE  if  the  MSEL 

model  is  to  be  used.  Otherwise,  it  is  FALSE. 

39.  KBAR/integer/ . Indicator  variable.  KBAR  =  1  if 

the  area  search  scenario  is  used.  KBAR  *  2  if  the 
barrier  search  scenario  is  used. 

40.  KCZO/integer/ . The  number  of  the  most  distant  CZ 

in  which  the  searcher  can  currently  make  a  detection. 


41.  KCZT/integer/ . The  number  of  the  most  distant  CZ 

in  which  the  target  can  currently  make  a  detection. 

If  KCZT  *  0,  no  CZ  detections  are  currently  possible. 
If  the  target  makes  a  CZ  detection,  KCZT  =  -1. 

42.  KP/*/integer/ . The  number  of  the  searcher  track 

anchor  point  at  which  the  search  pattern  begins  to 
repeat.  This  is  the  "return  point"  from  which 
subsequent  patterns  will  start  after  the  first  full 
pattern  is  executed  (input).  RESTR:  positive,  less 
than  NP.  See  NEXT (I)  and  NP. 

43.  MO/*/integer/ . The  number  of  range/propagation 

loss  data  points  used  to  define  the  searcher  direct- 
path  (DP)  propagation  loss  curve  (input) . 

RESTR:  greater  than  1,  and  less  than  21. 

44.  MODE/integer/ . Indicator  variable.  If  MODE  «  1, 

searcher  is  at  drift  speed.  If  MODE  *  2,  searcher 
is  at  sprint  speed. 

45.  MSELO/ integer/ . The  MSEL  counter  for  the  searcher. 

46.  MSELT/ integer/ . The  MSEL  counter  for  the  target. 

47.  MT/*/integer/ . The  number  of  range/propagation 

loss  data  points  used  to  define  the  target  direct 
path-propagation  loss  curve  (input) . 

RESTR:  greater  than  1,  and  less  than  21. 

48.  NBINS/*/integer/ . The  number  of  sections,  each 

with  a  fixed  probability,  that  define  the  initial 
target  distribution  across  the  top  of  a  barrier  choke 
point  (input) .  RESTR:  greater  than  1,  and  less 
than  21. 

49.  NBOTH/integer/ . The  number  of  simultaneous 

detections . 

50.  NCZDO (I) /integer/ . In  the  main  program,  the  number 

of  searcher  detections  made  in  the  Ith  convergence 
zone.  Converted  in  subroutine  SINKEM  to  the  fraction 
of  searcher  detections  in  the  Ith  CZ. 

51.  NCZDT( I) /integer/ . In  the  main  program,  the  number 

of  target  counter-detections  made  in  the  Ith  conver¬ 
gence  zone.  Converted  in  subroutine  SINKEM  to  the 
fraction  of  target  detections  in  the  Ith  CZ. 


52.  NCZO/*/integer/ . The  number  of  searcher  convergence 

zones  (input).  RESTR:  non-negative,  less  than  6. 

53.  NCZT/*/integer/ . The  number  of  target  convergence 

zones  (input).  RESTR:  non-negative,  less  than  6. 

54.  NDO/integer/ . The  number  of  detections  made  by 

searcher. 

55.  NDT/integer/ . The  number  of  detections  made  by  the 

target . 

56.  NEXT (I) /integer/ . The  next  anchor  point  in  the 

searcher  track  after  point  I.  The  program  sets 
NEXT (I)  *  1+1  for  I  *  1  to  NP-1,  and  NEXT(NP)  -  KP. 

Thus,  after  one  complete  execution  of  the  search 
pattern,  the  searcher  transits  from  point  NP  to  KP, 
and  repeats  the  pattern  starting  at  KP. 

57.  NL/ integer/ . Pointer  to  the  next  searcher  track 

leg  (similar  to  NLEG) . 

58.  NLEG/integer/ . A  pointer  to  indicate  the  searcher 

position  in  the  search  pattern.  The  searcher  is 
always  between  anchor  point  NLEG  and  point  NEXT (NLEG) . 

59.  NLS/*/integer/ . The  number  of  relative  bearing/ 

figure-of-merit  data  points  used  to  define  the  radial 
distribution  of  target  source  level  resulting  in  a  . 
radial  distribution  of  searcher  figure-of-merit 
(input).  RESTR:  greater  than  2,  less  than  51. 

60.  NONE/integer/ . The  number  of  replications  in  which 

neither  ship  makes  a  detection  prior  to  TMAX. 

61.  NP/*/integer/ . The  number  of  searcher  track  anchor 

points  (input).  RESTR:  greater  than  1,  less  than 
51. 

62.  NREP/*/interger/ . The  number  of  replications  to  be  run 

(input).  RESTR:  positive,  less  than  5001. 

63.  NTDO  Cl)  /real/hours/ . The  time  of  the  Ith  searcher 

detection . 

64.  NTDT (I) /real/hours/ . The  time  of  the  Ith  target 

counter-detection . 

65.  OFOMCI) /real/decibels/ . Through  the  EQUIVALENCE 

statement,  OFOMCI)  *  FOMOD,  0F0M(2)  =  FOMOS.  The 
subscript  of  this  variable  in  the  program  is  MODE. 

See  MODE,  FOMOD,  FOMOS. 


OL(I) /*/real/decibels/ . Direct-path  propagation 

loss  at  range  RO(I)  for  the  searcher  (input) . 
RESTR:  OL(I)  must  increase  as  range  increases. 
OL(I)  <  0LCI+1). 


PBIN  (I)  /*/real/ . The  probability  of  target  inclu¬ 

sion  in  the  Ith  section  when  the  user  defines  initial 
lateral  target  distribution,  I  =  1,...,NBINS,  (input). 
RESTR:  non-negative,  less  than  or  equal  to  1.0. 

Sum  of  all  PBIN (I)  =  1.0. 


RCZO(I, J) /*/real/nautical  miles/ . Range  to  the 

inner  (J=l)  and  outer  (J=2)  edges  of  the  Ith  conver¬ 
gence  zone  for  the  searcher,  1=1, . . . ,NCZO,  (input). 
RESTR:  RCZO(I,l)  <  RCZ0(I,2)  <  RCZO( 1+1,1)  < 

RCZO( 1+1,2) ,  etc. 


RCZT(I, J) /*/real/nautical  miles/ . Range  to  the 

inner  (J=l)  and  outer  (J=2)  edges  of  the  Ith  conver¬ 
gence  zone  for  the  target ,  I  =  1,...,NCZT,  (input). 
RESTR:  RCZT(I,1)  <  RCZT(I,2)  <  RCZT(I+1,1)  < 

RCZT(I+1,2),  etc. 


RFLCT/logical/ . RFLCT  is  FALSE  unless  a  target 

reflection  takes  place. 


RMAX/real/nautical  miles/ . The  maximum  of  RNGO, RNGT 


RNG/real/nautical  miles/ . Current  range  from 

searcher  to  target. 


RNGO/ real/nautical  miles/ . Searcher  current  direct 

path  detection  range. 


RNGT/real/nautical  miles/ . Target  current  direct 

path  detection  range. 


RNTDO (I ) /real/nautical  miles/ . Range  between 

searcher  and  target  when  searcher  detects  target. 


RNTDTd)  /real/nautical  miles/ . Range  between 

searcher  and  target  when  target  detects  searcher. 


R0(I) /*/real/nautical  miles/ . Range  associated 

with  searcher  propagation  loss  OL(I)  (input) . 
RESTR:  Must  be  entered  in  ascending  order,  RO(I)< 
R0CI+1)  . 


78.  RT(I)/*/real/nautical  miles/ . Range  associated 

with  target  propagation  loss  TL(I)  (input) . 

RESTR:  Must  be  entered  in  ascending  order,  RT(I)  < 
RT(I+1) . 

79.  RTCC/*/real/inverse  hours/ . Rate  of  target  course 

changes.  1/RTCC  is  the  mean  of  an  exponential  distri 
bution  from  which  the  time  interval  to  target  course 
change  is  drawn  (input) .  RESTR:  positive. 

80.  RTSC/*/real/inverse  hours/ . Rate  of  target  speed 

changes.  1/RTSC  is  the  mean  of  an  exponential  distri 
bution  from  which  the  time  interval  to  target  speed 
change  is  drawn  (input) .  RESTR:  positive. 

81.  SEED/*/integer/. ... .Psuedo- random  number  generator 
seed  (input) .  RESTR:  non-negative. 

82.  SIGMA (I) /*/real/decibels/ . Scale  parameter  for  the 

Ith  acoustic  fluctuation  process,  1=1, 2, 3. 

RESTR:  non-negative . 

83.  SINX/real/ . The  y-component  of  a  unit  vector  in 

the  direction  of  target  motion. 

84.  SO/real/knots/ . Current  searcher  speed. 

85.  SOD/*/real/knots/ . Searcher  drift  speed  (input) . 

For  constant  speed  search,  set  drift  speed  equal  to 
sprint  speed  (SOD  =  SOS) . 

86.  SOS/*/real/knots/ . Searcher  sprint  speed  (input) . 

For  constant  speed  search,  set  sprint  speed  equal  to 
drift  speed  (SOS  =  SOD)  . 

87.  SPECX/logical/ . This  variable  is  TRUE  if  the  user 

has  specified  a  lateral  target  distribution  other 
than  uniform  across  the  barrier  (barrier  scenario 
only).  Otherwise,  it  is  FALSE. 

88.  SPEED (I) /real/knots/ . Through  the  EQUIVALENCE 

statement,  SPEED (1)  =  SOD,  SPEED (2)  =  SOS.  The 
subscript  of  this  variable  in  the  program  is  MODE. 

See  MODE,  SOS,  SOD. 

89.  ST/real/knots/ . Target  current  speed,  drawn  from 

a  uniform  distribution  on  the  interval  (STMIN , STMAX) . 
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STACKl , STACK2 , STACK3 , STACK 4/real/ . Each  is  a  stack 

of  5000  psuedo-random  numbers  which  are  sequentially 
dereferenced  for  program  use  by  the  stack  counters 
C1,C2,C3,C4.  STACKl  is  distributed  exponentially, 
mean  *  1.0.  STACK2  is  sistributed  normally,  mean  - 
0.0,  variance  *  1.0.  STACK3  is  distributed  uniform- 
ally  on  the  interval  (0,1).  STACK4  is  distributed 
as  specified  by  the  user. 


START/*/real/hours/ . The  latest  time-late  on 

barrier  time.  The  time-late  on  barrier  for  the 
searcher  is  distributed  uniformly  on  the  interval 
(0, START),  with  an  assumption  that  the  target  travels 
at  mean  speed  from  time  zero  until  the  searcher  is  on 
the  barrier  (input).  RESTR:  non-negative. 


STATSl/logical/ . This  variable  is  FALSE  if  no 

statistical  analysis  is  to  be  done.  Otherwise  it 
is  TRUE. 


STATS2/logical/ . This  variable  is  TRUE  if  section¬ 

ing  of  the  data  is  to  be  accomplished.  Otherwise, 
it  is  FALSE. 


STINC/real/knots/ . The  interval  of  speed  which 

defines  the  uniform  distribution  of  allowable  target 
speeds  (STMAX  -  STMIN) . 


95.  STMAX/*/realAnots/ . Target  maximum  speed  (input)  . 

RESTR:  STMAX  must  be  greater  than  or  equal  to  STMIN. 


STMIN/*/realAnots/ . Target  minimum  speed  (input)  . 

RESTR:  STMIN  must  be  less  them  or  equal  to  STMAX. 


TBIG/real/nautical  miles  or  hours/ . TBIG  is  used 

to  compute  the  simulation  time  step.  Initially,  TBIG 
is  a  range  ■  RNG  -  RMAX,  and  is  positive  if,  and  only 
if,  no  direct  path  detection  is  possible.  If  no 
detection  is  possible  (CZ  or  DP) ,  TBIG  is  set  to  the 
minimum  closing  distance  such  that  either  a  detection 
or  counter-detection  is  just  possible  (DP  or  CZ) 
based  on  the  current  acoustic  conditions.  TBIG  is 
converted  to  a  time  interval  by  dividing  by  the 
current  combined  speeds  of  the  searcher  and  target, 
which  assumes  a  worst  case  head-to-head  closing 
situation.  The  simulation  time-step  is  then  set  to 
the  maximum  of  TBIG  or  0.05  hour. 


98.  TCC/realAours/ . Time  of  next  searcher  course 

change . 


99.  TD/*/real/hours/. 
RESTR:  positive . 


Searcher  drift  time  (input) . 


100.  TPOM (I) /real/decibels/ . Through  the  EQUIVALENCE 

statement,  TFOM(l)  =  FOMTD,  TFOM(2)  *  FOMTS.  The 
subscript  of  this  variable  in  the  program  is  MODE. 

See  MODE,  FOMTD,  FOMTS. 

101.  THETA/real/radians/ . Target  psuedo-course .  THETA 

is  measured  from  the  x-axis  in  a  counter-clockwise 
direction. 

102.  TIFL (I) /real/hours/ . For  1=1, 2, 3,  TIFL(I)  is  the 

time  at  which  AFL(I)  changes  in  the  LSJ  process,  or 
the  next  sample  time  for  the  GMA  process.  For 

I  *  4-7,  through  the  EQUIVALENCE  statement, 

TIFL (4)  =  TCC,  TIFL (5)  =  TTCC,  TIFL (6)  =  TTSC, 

TIFL (7)  =  TSC. 

103.  TIME (I) /real/hours/ . Through  the  EQUIVALENCE 

statement,  TIME (11  =  TD,  TIME (2)  =  TS.  The  subscript 
of  this  variable  in  the  program  is  MODE.  See  MODE, 
TD,  TS. 

104.  TINC/real/hours/ . The  time  at  which  the  next 

simulation  time-step  ends. 

105.  TL (I) /*/real/decibels/ . Propagation  loss  at  range 

RT (I)  for  the  target  (input).  TL(I)  must  increase 
with  increasing  range.  TL(I)  must  be  less  than 

TL  (I  +  1)  . 

106.  TLAST/ real/hours/ . The  time  at  which  the 

Gauss-Markov  error  process  was  last  evaluated. 

107.  TMAX/*/real/hours/ . Maximum  allowed  search  time 

per  replication. 


108. 

109. 


110. 


111. 


112. 


TNOW/real/hours/ . Current  time  in  a  replication. 


TO  (I) /real/hours/ . The  time  at  which  the  figure- 

of-merit  sample  for  the  searcher,  FO(I) ,  was  taken. 


TOXMAX/real/nautical  miles/ . Twice  the  value  of 

XMAX. 


TOYMAX/ real/nautical  miles/ . Twice  the  value  of 

YMAX. 


TS/*/real/hours/ . Searcher  sprint  time  (input) . 

RESTR:  positive. 
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TSC/real/hours/ . Time  of  next  searcher  speed 

change . 

TT( I) /real/hours/ . The  time  at  which  the  target 

figure-of-merit  sample,  FT(I) ,  was  taken. 

TTCC/real/hours/ . Time  of  next  target  course 

change. 

TTSC/real/hours/. . . . .Time  of  next  target  speed 
change . 

TWOPI/real/ . Twice  the  value  of  it. 

UX/real/knots/ . The  x-component  of  current  target 

velocity  vector. 

UY/real/knots/ . The  y-component  of  current  target 

velocity  vector. 

VX/ real/knots/ . The  x-component  of  current 

searcher  velocity  vector. 

VY/rea 1/knots/ . The  y-component  of  current 

searcher  velocity  vector. 

XBIN  (I )  / Vreal/nautical  miles/ . The  distance  from 

the  origin  (x  *  0.0)  to  the  right-most  x-value  of 
the  section  containing  target  probability  mass 
PBIN(I).  RESTR:  positive,  XBIN(NBINS)  =  XMAX. 

XMAX/*/real/nautical  miles/ . Length  of  the  search 

area  rectangle  in  the  x-direction  (input) . 

RESTR:  positive . 

XO/real/nautical  miles/ . The  current  searcher 

x-position. 

XODT CD /real/nautical  miles/ . The  x-position  of 

the  searcher  when  the  searcher  makes  the  Ith 
detection  of  the  target. 

XOT/rea 1/nautical  miles/ . The  projected  searcher 

position  based  on  the  current  position  (XO)  and  the 
computed  simulation  time-step. 

XOTD (I) /real/nautical  miles/ . The  x-position  of 

the  searcher  when  the  target  makes  the  Ith  counter¬ 
detection  of  the  searcher. 


128.  XP (I) /*/real/nautical  miles/ . The  x-position  of 

the  Ith  searcher  track  anchor  point  (input) . 

RESTR:  non-negative. 

129.  XT/real/nautical  miles/ . The. current  x-position 

of  the  target. 

130.  XTDO( I) /real/nautical  miles/ . The  x-position  of 

the  target  when  the  target  makes  the  Ith  counter¬ 
detection  of  the  searcher. 

131.  XTOD(I) /real/nautical  miles/ . The  x-position  of 

the  target  when  the  searcher  makes  the  Ith  detection 
of  the  target. 

132.  YMAX/*/real/nautical  miles/ . Length  of  the  search 

area  in  the  y-direction  (input)) . 

RESTR:  positive. 

133.  YO/real/nautical  miles/ . Current  searcher 

y-position. 

134.  YODT (I) /real/nautical  miles/ . The  y-position  of 

the  searcher  when  the  searcher  makes  the  Ith 
detection  of  the  target. 

135.  YOT/real/nautical  miles/ . Projected  searcher 

y-position  based  on  the  current  position  (YO)  and 
the  simulation  time-step. 

136.  YOTD( I) /real/nautical  miles/ . The  y-position  of 

the  searcher  when  the  target  makes  the  Ith  counter¬ 
detection  of  the  searcher. 

137.  YP (I) /*/real/nautical  miles/ . The  y-position  of 

the  Ith  searcher  track  anchor  point  (input) . 

RESTR:  positive . 

138.  YT/real/nautical  miles/ . Current  y-position  of 

the  target. 

139.  YTDO (I) /real/nautical  miles/ . The  y-position  of 

the  target  when  the  target  makes  the  Ith  counter¬ 
detection  of  the  searcher. 

140.  YTDO (I) /real/nautical  miles/ . The  y-position  of 

the  target  when  the  searcher  makes  the  Ith  detection 
of  the  target. 


ABBREVIATIONS  AND  ACRONYMS 


This  appendix  contains  an  alphabetical  listing  of  abbrev 
iations  and  acronyms  used  in  this  thesis.  . 

1.  ASW:  Antisubmarine  Warfare. 

2.  Co v ( X ,  Y ) :  If  X, Y  are  random  variables,  this  is  the 
covariance  function.  If  X,Y  are  realizations  of  an 
autoregressive  time  series,  this  is  the  autocovari¬ 
ance  function. 

3.  CPU:  Central  Processing  Unit. 

4.  CZ:  Convergence  Zone.  Usually  used  to  describe  the 
propagation  mode  of  an  acoustic  signal. 

5.  db:  Decibel. 

6.  DP:  Direct  Path.  Usually  used  to  describe  the  propa¬ 
gation  mode  of  an  acoustic  signal. 

7.  E (X) :  The  expected  value  of  the  random  variable  X. 

8.  ER:  Exchange  Ratio. 

9 .  EXP (x) :  An  exponential  distribution  with  rate 
parameter  x. 

10.  FOM:  Figure-of-Merit. 

11.  IMSL:  International  Mathematics  and  Statistics 
Library . 

12.  Le:  Environmental  noise. 

13.  Ls:  Target  noise  source  level. 

14.  MOEs:  Measure  Of  Effectiveness  for  the  Searcher 

15.  MOEt:  Measure  Of  Effectiveness  for  the  Target. 

16.  MSEL:  Minimum  Signal  Excess  Logic. 


N(x,y) :  A  normal  distribution  with  mean  x  and 
variance  y. 

NON IMS L:  Non- International  Mathematics  and  Statistics 
Library . 

NPGS(NPS):  Naval  Postgraduate  School. 

Nrd:  Recognition  differential. 

PASS:  Passive  Acoustic  Search  Simulation. 

PCD:  Probability  of  Counter-Detection. 

PCDCZ :  Probability  of  Counter-Detection  in  a 
Convergence  Zone . 

PCDCZI:  Probability  of  Counter- Detect ion  in  the  Ith 
Convergence  Zone. 

PCDDP:  Probability  of  Counter-Detection  by  Direct 
Path  propagation. 

PD:  Probability  of  Detection. 

PDCZ:  Probability  of  Detection  in  a  Convergence  Zone. 

PDCZI:  Probability  of  Detection  in  the  Ith 
Convergence  Zone. 

PDDP:  Probability  of  Detection  by  Direct  Path 
propagation . 

p(X,Y) :  If  X,Y  are  random  variables,  p  is  the 
linear  correlation  coefficient.  If  X,Y  are 
realizations  of  an  autoregressive  time  series, 
then  p  is  the  autocorrelation  function. 

Rs/Rt:  The  mean  detection/counter-detection  range. 

(R50) s/(R50) t:  The  median  detection/counter¬ 
detection  range. 

SE:  Signal  Excess. 

Ts/Tt:  Mean  time  to  detection/counter-detection. 

(T50) s/CT50) t:  Median  time  to  detection/counter¬ 
detection. 
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APPENDIX  D 

APPLICATION  OF  THE  PASSIVE  SONAR  EQUATION  TO  THE  PASS  MODEL 


The  passive  sonar  equation  may  be  written  as : 


SNR  =  Ls  -  Nw  -  Le 


where : 


1.  SNR  *  the  signal-to-noice  ration  at  the  processor 
output  (db) . 

2.  Ls  *  the  target  source  level  (db) . 

3.  Nw  =  the  signal  propagation  loss  (db) . 

4.  Le  *  the  background  noise  compensated  for  by  sonar 
directivity  (db) . 

Now,  make  the  following  definitions: 

1.  Nrd  *  RECOGNITION  DIFFERENTIAL  (db) .  Nrd  is  that  SNR 
required  for  a  probability  of  detection  of  0.5.  It 
can  be  thought  of  as  the  SNR  required  for  a  target  to 
be  recognized  as  a  valid  contact  50%  of  the  time.  Nrd 
is  generally  a  function  of  sonar  processor  performance 

2.  FOM  *  FIGURE-OF-MERIT  (db) .  FOM  is  that  amount  of 
signal  propagation  loss  (Nw)  that  results  in  a 
probability  of  detection  of  0.5. 

3.  SE  =  SIGNAL  EXCESS  (db) .  SE  is  the  algebraic  differ¬ 
ence  of  FOM  and  Nw  (SE  *  FOM  -  Nw) . 

If  we  now  substitute  Nrd  for  SNR  and  FOM  for  Nw  in  the  above 

equation,  and  solve  for  FOM,  we  have: 


FOM  ■  Ls  -  Le  -  Nrd 


Thus,  the  result  is: 


SE  >  0 
SE  <  0 


Nw  <  FOM 
Nw  >  FOM 


Pd  >  0.5 
Pd  <  0.5 


In  PASS,  we  obtain  a  certain  detection  whenever  the  sum 
of  the  deterministic  signal  excess  and  the  error  term  is 
greater  than  or  equal  to  zero.  TOie  error  term  has  a  zero 
mean,  and  the  probability  that  it  is  greater  than  zero  is 
0.5.  Thus,  when  the  deterministic  signal  excess  is  zero,  the 
probability  of  detection  in  PASS  is  0.5,  which  is  consistent 
with  the  passive  sonar  equation  results  above. 


APPENDIX  E 


COMPUTER  SIMULATION  OF  THE  STOCHASTIC  ERROR  PROCESSES 

In  this  appendix  we  examine  the  modeling  of  the  acoustic 
fluctuation  process  independent  of  all  other  events  in  the 
simulation.  In  effect,  we  will  ignore  all  events  which 
affect  the  simulation  time-step  except  for  acoustic 
fluctuations . 

A.  THE  LAMBDA-SIGMA  JUMP  PROCESS 

In  PASS  there  are  three  fluctuation  processes  simultane¬ 
ously  in  progess. 

X(t)  ■  searcher  local  error  process. 

Y (t)  ■  global  error  process 

Z(t)  *  target  local  error  process 

Each  process  is  defined  by  a  rate  parameter  (X^,  X2 ,  X3) 
and  a  scale  parameter  (o^,  a 2 ,  a3).  For  illustrative 
purposes,  assume  we  initially  start  the  process  at  time  *  T 
with  independent  draws  from  normal  distributions  parameter¬ 
ized  by  the  scale  parameters.  That  is: 

,  2. 

XtT)  »  where  nx  ~  N  (o,^) 

Y(t)  ■  riy  where  ny  ~  N  <o,a2) 

Z(t)  *  n2  where  nz  ^  N  (o,a|) 
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Also,  at  time  *  T,  we  determine  the  time  interval  to  the 
next  fluctuation  level  change,  s,  by  independent  draws  from 
exponential  distributions  parameterized  by  the  rate 
parameters .  That  is : 


Sx  "  Tx  where  Tx  *  EXP ( A L) 


Sy  =  ty  where  ty  ^  EXPU2> 


Sz  =  tz  where  t2  'v.  EXP(X3) 


The  actual  time  of  the  fluctuation  changes  would  then  be: 


t  *  S  +  T 
xx 


t  *=  S  +  T 

y  y 


t  =  S  +  T 
z  z 


If  no  detection  takes  place  at  time  =  T,  the  next  time  a 
detection  can  take  place  is  the  time  of  the  first  fluctuation 
change.  Therefore,  the  simulation  time-step  ,  At,  would  be: 

At  =  min  {Sx,Sy,Sz> 

All  simulation  events  (e.g.  platform  motion)  would  then  be 
executed  with  time-step  At,  and  the  simulation  time  would  be 
incremented  by  At. 

The  fluctuation  level  associated  with  At  is  then 
changed.  Suppose,  for  example,  that  tx<ty<tz*  Then,  we 


would  change  the  X  fluctuation  level  by  replacing  the 

2 

previous  level  with  a  draw  from  N(0,a^),  and  obtain  a  new 
time  interval  to  the  next  change  of  X  by  making  a  draw  from 
EXP(A^) ,  from  which  the  time  of  the  next  change  in  X  is 
calculated.  If  no  detection  takes  place,  a  new  time-step  is 
computed  as  before,  and  the  process  is  repeated. 

B.  THE  GAUSS -MARKOV  PROCESS 

Unlike  the  LSJ  process,  the  GMA  process  involves  continu¬ 
ous  sample  paths.  The  defining  function  of  the  GMA  process 
is  of  the  form: 

Xtt  +  s)  «  e“Xs  X(t)  +  g(s)nx 

where : 

1.  n  is  a  normally  distributed  random  variable  with 

2 

zero  mean  and  variance  o^. 

2.  g(s)  is  a  function  of  the  time  increment,  s,  which 
has  the  properties: 

lim  g(s)  *  0  and  lim  g(s)  =  1 
s-*o  s-*<» 

we  know  that: 

Var [X(t) ]  -  a*  Vt 

so 

Var  [X(t+s>  ]  -  a\  -  Var  [e"Xs  X(t)  +  g(s)nx3 


Since  e  and  g(s)  are  not  random  variables  and  X(t)  and  n: 
are  independent,  we  can  write: 

a 2  -  e’2Xs  Var[X(t)]  +  CgCs)]2  Var[nxJ 


a2  =  e-2Xs  g2  +  [g(g)]2  a2 


g(s)  =  (l-e“2Xs)1/2 


Thus ,  the  functional  form  of  GMA  is: 

XCt+s)  *  pXCt)  +  U-P2)1/2  nx 
where: 


The  difficulty  in  using  this  process  in  a  computer  simu¬ 
lation  is  deciding  on  the  time-step.  The  more  accurately  one 
wants  to  model  the  process,  the  smaller  the  time-step 
required,  but  small  time-steps  cause  excessively  long  run 
times.  As  a  practical  compromise,  the  PASS  model  uses  the 
identical  method  for  determining  fluctuation  time-steps  in 
both  the  LSJ  and  GMA  process  (i.e.  the  GMA  process  is 
sampled  at  exponentially  distributed  random  times  governed 
by  the  rate  parameters) .  The  difference  in  applying  this  to 
the  GMA  process  is  that  each  time  one  fluctuation  process  is 
changed  in  LSJ,  all  three  are  updated  in  GMA. 


For  example,  in  the  previous  section  where  X  was  changed, 
the  GMA  process  would  update  X,  Y,  and  Z  based  on  a  time 


increment  s  ,  and  independent  draws  to  a,  ti  /  n,« 
x  x  y  « 

C.  GENERATION  OF  PSUEDO  RANDOM  NUMBERS 

Psuedo  random  numbers  used  throughout  PASS  are  generated 
using  the  NPGS  random  number  package  LLRANDOMII .  This 
package  was  used,  vice  a  portable  package  (e.g.  International 
Mathematics  and  Statistics  Library)  because  o£  the  improve¬ 
ment  in  speed  in  number  generation.  The  IBM  Assembly 
Language  LLRANDOMII  routines  are  between  three  and  eight 
tiroes  as  fast  as  the  FORTRAN  IMSL  subroutines. 

D.  CORRELATION  OF  THE  ERROR  PROCESSES 

Suppose  we  define  two  compound  error  processes,  A  and  B, 
as  follows: 

Att)  =  XCt)  +  Ytt) 

BCt)  *  Ytt)  +  ZCt) 

where  X,  Y,  Z  are  LSJ  or  GMA  processes.  Dropping  the  ”t" 
notation  for  convenience,  a  computational  formula  for  the 
covariance  is: 

Cov(A,B)  *  E[A*B]-EIA]  *E[B] 

and  since  E[A]  *  E[B]  >0 

Cov(A,B)  *  E [A*B]  -  E[XY+Y2+XZ+YZ] 
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Cov(A,B)  *  Cov(X,Y)  +  a2  +  Cov (X, Z)  +  Cov(Y,Z) 
Cov(A,B)  *  a2 

since  X,Y,Z  are  independent,  and  E[Y]  =  0. 

The  linear  coefficient  of  correlation  is  defined  as: 

p(A,B>  =  522%SL 


where 


B 


Co*  +  a*)1'2 

,2  ,  2.1/2 

(o2  +  a3) 


so,  the  coefficient  of  correlation  for  the  compound  error 
process  is: 


Note  that:  1.  lim  ptA,B)  *  1 


which  satisfies  the  intuitive  and  practical  requirement  that 
as  the  global  fluctuation  becomes  dominant,  the  error  signal 
experienced  by  both  sensors  will  be  highly  correlated. 


2.  lim  p(A,B)  *  0 


which  satisfies  the  intuitive  and  practical  requirement  that 
as  the  local  fluctuations  dominate,  the  error  signal  experi¬ 
enced  by  the  sensors  will  be  almost  uncorrelated. 
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A  THREE-OUT-OF-FIVE  DETECTION  CRITERIA  MODEL 

The  Minimum  Signal  Excess  Logic  (MSEL)  starts  in  state 
zero  (MSEL  0) .  Each  time  the  acoustic  environment  is 
samples,  the  MSEL  state  is  incremented  by  one  if  detection 
is  possible  (SE  >  0) .  If  detection  is  not  possible  (SE  <  0) , 
then  the  MSEL  counter  is  decremented  by  one,  unless  MSEL  is 
in  state  zero.  The  lowest  MSEL  state  is  zero,  and  detection 
occurs  the  first  time  the  MSEL  state  reaches  three.  This 
logic  results  in  a  detection  occur ing  if,  and  only  if,  three 
of  the  last  five  samples  of  :the  acoustic  environment  had 
SE  2  0. 

MSEL  is  equivalent  to  the  Markov  chain  shown  in  Figure 
F.l,  where  the  transition  probabilities  are,  in  general,  a 
function  of  time  (i.e.  they  will  change  each  time  the 
acoustic  signal  is  sampled) .  in  using  the  LS J  process ,  this 
presents  no  problem  in  that  the  fluctuation  process  is 
static  between  samples.  In  using  the  GMA  process,  however, 
the  fluctuation  process  is  dynamic  between  sample  times,  and 
therefore  the  rate  of  sampling  will  have  an  effect  on  the 
time  of  detection.  The  PASS  model  accounts  for  this  in  two 
ways ; 

1.  The  rate  of  environment  sampling  is  the  same  for  LSJ 
and  GMA.  See  Appendix  E. 

2.  Whenever  either  searcher  or  target  are  not  in  MSEL  0, 
the  maximum  PASS  time-step  defaults  to  0.05  hour  (3 
minutes) . 


For  analysis  purposes,  assume  all  the  transition  prob¬ 
abilities  are  constant  and  equal  to  p.  Let  q  =  (1-p)  and 
let  Si  =  the  mean  number  of  time  increments  (samples  of  the 
acoustic  environment)  to  go  from  state  i  to  state  3.  Reach¬ 
ing  the  absorbing  state,  state  3,  is  equivalent  to  detection. 

We  can  write  a  series  of  balance  equations  as: 

S2  =  p(l)  +  q(l  +  SL)  =  1  +  qS^ 

=  p(l  +  S2)  +  q (1  +  So)  *  1  +  P^2  +  qSo 

So  —  p(l  +  Sx)  +  q (1  +  So)  *  1  +  pSL  +  qSo 

If  we  solve  this  system  of  equations  for  So,  we  have: 

So  =  |(2  +  iy) 

F  P 

Since  So  is  the  mean  number  of  time  steps  to  go  from 
MSEL  state  zero  to  MSEL  state  3  (for  a  fixed  p) ,  we  can 

place  an  upper  bound  on  the  expected  time  to  go  from  MSEL  0 

to  MSEL  3  as: 

To  =  3So  (in  minutes) 

since  the  maximum  time-step  is  PASS  is  three  minutes  when 
the  searcher  or  target  MSEL  counter  is  greater  than  zero. 
Table  VIII  shows  representative  values  of  So  and  To  as  a 
function  of  the  instantaneous  probability  of  detection. 
Keeping  in  mind  the  simplifying  assumptions  made,  the  values 
of  To  are  analogous  to  maximum  expected  integration  times. 


a 


Figure  F. 

1.  The  MSEL  Model  as  a  Markov  Chain. 

TABLE  VIII 

MSEL 

MAXIMUM  EXPECTED  INTEGRATION 

TIMES 

P 

SO 

To 

1.0 

3.000 

9.00 

0.9 

3.594 

10.78 

o 

. 

00 

4.453 

13.36 

0.7 

5.773 

17.32 

VO 

• 

o 

7.963 

23.89 

in 

• 

o 

12.000 

36.00 

0.4 

20.625 

61.88 

0.3 

43.704 

131.11 

0.2 

135.000 

405.00 

0.1 

1020.000 

3060.00 

iWi 
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CODE  FOR  EXECUTIVE  PROGRAM  TO  RUN  PASS 


&TRACE 

CP  TERMINAL  LINES I ZE  80 

ERASE  PASS  OUTPUT 

ERASE  PASS  HOLD 

ERASE  PASPLT  DATA 

COPYFILE  PASS  DATA  Al  =  HOLD  = 

ERASE  PASS  DATA 

FILEDEF  06  DISK  PASS  OUTPUT 

FILEDEF  07  DISK  PASS  DATA 

FILEDEF  08  TERM 

FILEDEF  09  DISK  PASS  HOLD 

FILEDEF  10  DISK  PASPLT  DATA 

CLRSCRN 

STYPE  PASSIVE  ACOUSTIC  SEARCH  SIMULATION  ***  "PASS" 

&  TYPE 

STYPE  EXECUTING  PROGRAM  "PASS"  **  "PASS"  MUST  BE  FORTRAN-G 
OR  FORTRAN-H 

STYPE  PRECOMPILED  ON  YOUR  A-DISK.  THE  NONIMSL  LIBRARY  MUST 
BE  AVAILABLE. 

STYPE 

STYPE  DATA  FILE  MANAGEMENT: 

&TYPE 

STYPE  A.  INPUT  DATA  FILE  "PASS  DATA": 

STYPE  THIS  FILE  CONTAINS  THE  MOST  RECENT  DATA  FOR  RUNNING  PASS. 
STY RE  IT  IS  CREATED,  OR  MODIFIED,  AFTER  RUNNING  THE  INTERACTIVE 
STYPE  DATA  INPUT  SECTION  OF  PASS.  SEE  FURTHER  NOTES  ON  OUTPUT 
DATA 

STYPE  FILE  "PASS  HOLD". 

STYPE 

STYPE  B.  OUTPUT  DATA  FILE  "PASS  HOLD": 

STYPE  UPON  EXECUTIO  OF  PASS,  IF  FILE  "PASS  DATA"  ALREADY 
EXISTS 

STYPE  IT# IS  COPIED  INTO  "PASS  HOLD"  UNCHANGED.  THUS  IF  YOU 
CHANGE  THE 

STYPE  INPUT  DATA  DURING  THE  INTERACTIVE  SECTION,  THE  CHANGED 
DATA  WILL 

STYPE  BE  IN  "PASS  DATA",  AND  THE  UNCHANGED  DATA  WILL  BE  IN 
"PASS  HOLD". 

STYPE 

STYPE  C.  OUTPUT  DATA  FILE  "PASPLT  DATA": 

STYPE  THIS  FILE  CONTAINS  ALL  THE  DATA,  IN  THE  PROPER  FORMAT, 
NEEDED 

STYPE  TO  GENERATE  GRAPHICAL  OUTPUT  USING  PROGRAM  "PASSPLT" 

AND  THE 
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STYPE  DISSPLA  GRAPHICS  SYSTEM. 

&TYPE 

&TYPE  D.  OUTPUT  DATA  PILE  "PASS  OUTPUT": 

6TYPE  THIS  PILE  CONTAINS  THE  STATISTICAL  OUTPUT  FROM  THE 
MOST  RECENT 
&TYPE  RUN  OF  PASS. 

STYPE 

6 TYPE  - 

LOAD  PASS  (START 
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FORTRAN  CODE  FOR  PASS 
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